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We examine n component spin systems with arbitrary two spin interactions (of unspecified range) 
within a general framework to highlight some new subtleties present in incommensurate systems. We 
determine the ground states of all translationally invariant 0(n > 2) systems and prove that barring 
commensurability effects they are always spiral- no other ground states are possible. We study the 
effect of thermal fluctuations on the ground states to discover a novel odd-even n effect. Soft spin 
analysis suggests that algebraic long range order is possible in certain frustrated incommensurate 
even n systems while their odd n counterparts exhibit an exponential decay of correlations. We 
illustrate that many frustrated incommensurate continuous spin systems display smectic like ther- 
modynamics. We report on a generalized Mermin-Wagner-Coleman theorem for all two dimensional 
systems (of arbitrary range) with analytic kernels in momentum space. A new relation between gen- 
eralization Mermin-Wagner-Coleman bounds and dynamics is further reported. We suggest a link 
between a generalized Mermin-Wagner-Coleman theorem to divergent decoherence (or bandwidth) 
time scale in the quantum context. A generalization of the Peierls bound for commensurate systems 
with long range interactions is also discussed. We conclude with a discussion of 0(n) spin dynamics 
in the general case. 
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I. INTRODUCTION 

In this article (constituting a portion of my thesis [1]), 
we aim to unveil some of general properties of 0(n) spin 
systems having two-spin interactions. We investigate 
Ising and 0(n) ground states, derive generalized Mermin- 
Wagner inequalities and illustrate how Peierls' bounds 
may be derived in some long range systems. Unlike the 
very rich ground state structures possible for incommen- 
surate scalar systems (bubbles, dots, Wigner crystals), 
we prove here that for incommensurate continuous spin 
systems on a lattice, spiral states are the only ground 
states. Perhaps most noteworthy, we perform a thermal 
stability analysis whose results coincide with a general- 
ized Mermin- Wagner inequality. In the general case, we 
discover a new intriguing odd-even n effect that is ab- 
sent in conventional ferromagnetic or antiferromagnetic 
systems: Algebraic long range order may be possible for 
certain frustrated incommensurate even n systems, while 
their odd n counterparts exhibit an exponential decay of 
correlations. Odd n incommensurate systems are gener- 
ally more disordered than their even n counterparts. In 
particular, we illustrate that many frustrated incommen- 
surate continuous spin systems have smectic like thermo- 
dynamics. 

The outline of this article is as follows: After introduc- 
ing, in section(II), the terms and notations that will be 
used throughout the work, we provide, in section (III) 
frustrated toy models which are rich enough to illus- 
trate some of the general features that we aim to high- 
light. These models will be employed for illustrative 
purposes only. The contents of the present publication 



are not limited to these systems. All stated in the cur- 
rent article applies to all translationally invariant sys- 
tems (whether these are the standard nearest neighbor 
Heisenberg model or the more complex frustrated mod- 
els that we introduce in section (III)). Notwithstanding 
their pedagogical purposes in the current context, the 
toy models presented in section(III) and myriad variants 
have been a subject of much research in recent years, 
see e.g. [2]. The physics of incommensurate orders along 
with, in some instances, competing short- and long-range 
interactions is a rich and complex topic. These appear, 
amongst others, in the high temperature superconduc- 
tors [3], [4], [5], [6], [7], manganates and nickelates [8], 
[9], quantum Hall systems [10], [11], [12] chemical and 
magnetic mixtures, crumpled membranes [2], and some 
theories of structural glasses [13], [14]. In many instances, 
there appear to be a large number of possible ground 
states exhibiting nanoscale phase separation, including 
stripes, labyrinths, patches, and dot arrays. Most new 
results reported in this article pertain to incommensu- 
rate spin structures. 

In section(IV) we discuss the ground states of Ising 
spin models and show what patterns one should expect 
in general. Once the ground states will be touched on, we 
will head on to show how Peierls bounds may be estab- 
lished for many systems having infinite range interactions 
if the ground states are simple. In section(VI) we shortly 
review mean field solutions of the general two spin Ising 
models. A generalization of the standard Peierls bound 
to other systems with long range interactions is provided 
in section (V). When fused with an additional sym- 
metry (which is not present in many of the models that 
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we discuss), Peierls' bounds suffice to prove the existence 
of long range order. 

Henceforth, the bulk of this article focuses on con- 
tinuous spins whose number of components n > 2. In 
section (VII) we prove that, sans special commensura- 
bility effects, the ground states of all 0(n > 2) will 
typically have a spiral like structure. Section(VIII) de- 
tails an exceedingly simple spin wave stiffness analysis to 
gauge the effect of thermal fluctuations on the various 
0(n > 2) ground states. In section(IX) we will discuss 
thermal fluctuations within the framework of "soft-spin" 
XY model. We will see that the normalization constraint 
gives a Dirac like equation. In the aftermath, the fluctu- 
ation spectrum will be seen to match with that derived 
in section(VIII). We will show possible links to smectic 
like behavior in three dimensions. Next, we go one step 
further to study the fully constrained "hard-spin" 0(2) 
and 0(3) models and show (in section X) that all trans- 
lationally invariant systems in two dimensions with an 
analytic rotationally symmetric interaction kernel never 
develop spontaneous magnetization. At the end of the 
section our analysis will match that of sections(VIII) and 
(IX) We extend the Mermin-Wagner-Coleman bounds, in 
section(XI), to high dimensions to show that intricacies 
occur if a certain high dimensional integral will be seen 
to diverge. In the low temperature limit, the integrand of 
this integral will, once again, match that derived by the 
much more naive spin stiffness and soft spin fluctuation 
analysis of earlier sections. In section(XII) we will exam- 
ine the "soft-spin" version of Heisenberg spins. We will 
see that it might be naively expected that the spin fluc- 
tuations in odd n spin systems are larger than in those 
with an even number of spin components. The origin of 
this "odd" even-odd effect is that for odd n systems, one 
of the spin components is unpaired and may exhibit less 
inhibited fluctuations. Next, in section(XIII), we carry 
out the spin fluctuation analysis for four component soft 
spins to see that their spectra coincides with that pre- 
dicted in the earlier spin stiffness analysis. We compute 
the correlation functions for this system and find that 
even n incommensurate systems may display algebraic 
long range order in situations where their odd n coun- 
terparts of the same system exhibit exponential decay of 
correlations at all finite temperatures (insofar as we may 
discern from a perturbative analysis of soft spin models). 

We match the finite n analysis with its large n counter- 
part in section(XIV). We show that in the limit of large 
n both odd and even component spin systems behave in 
the same manner. Essentially, they all tend towards an 
"odd" behavior. We also report on a "holographic" like 
effect present in some frustrated systems. In these, the 
ground state entropy is shown to scale with the surface 
area of the system. In section(XVI) , we compute the crit- 
ical temperature of all translationally invariant 0(n > 2) 
spin models within mean field theory. In section(XVII) 
we briefly remark that much of analysis is unchanged 



for arbitrary non-translationally invariant two spin in- 
teractions. We conclude with a discussion of 0(n) spin 
dynamics. 

A central theme which will be repeatedly touched on 
throughout the paper is the possibility of non-trivial 
ground state manifolds. If the system is degenerate the 
effective topology of the low temperature phase of the 
system may be classified in momentum (or other basis). 
In such instances the low temperature behavior of the 
systems will be exceedingly rich. 

II. DEFINITIONS 

We consider simple classical spin models of the type 
H=lY.V&y)lS(x)-S(y)]. (1) 

Here, the sites x and y lie on a (generally hypercu- 
bic) lattice of size N. The spins {S(x)} are normalized 
and have n components, Sf(x) ~ 1> & t au lattice 

sites x. We will primarily focus on translationally in- 
variant interactions V(x, y) = V(x — y). We employ 
the non-symmetrical Fourier basis convention (f(k) — 
Y<x F (x) e ~ i% ' S i F W = wEfe/W^) wherein the 

Hamiltonian is diagonal and reads 

ff = ^E^)i 5? wi 2 ( 2 ) 

k 

where v(k) and S(k) are the Fourier transforms of V(x) 
and S(x). More generally, for some of the properties 
that we will illustrate, one could consider any arbitrary 
real two spin interactions (x|y|y) which would be diago- 
nalizcd in another basis {\u}} instead of the Fourier ba- 
sis. For simplicity, we will set the lattice constant to 
unity- i.e. on a hypercubic lattice (of side L) with pe- 
riodic boundary conditions the wave-vector components 
h = where n is an integer (and the real space co- 
ordinates xi are integers). Throughout, we employ K to 
denote reciprocal lattice vectors and A(fc) as a shorthand 
for the lattice lattice Laplacian: 

d 

A(jfe) = 2^(1 -cos jfej). (3) 
i=i 

In some of the frustrated systems that we will soon con- 
sider, v(k) may be written explicitly as the sum of sev- 
eral terms: those favoring homogeneous states k — > 0, 
and those favoring zero wavelength k — > oo (or k — > 
(tt, tt, 7r) on a lattice.) As a result of this competition, 
modulated structures arise on an intermediate scale. 
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III. TOY MODELS 

Although we will keep the discussion very general, it 
might be useful to have a few explicit applications in 
mind. There is a lot of physical intuition which under- 
lies the upcoming models. Unfortunately, here they will 
merely serve as nontrivial toy models on which we will 
able to exercise our newly gained intuition. The systems 
to be presented are frustrated: not all two spin interac- 
tions can be simultaneously satisfied. They entail com- 
peting interactions. Such systems appear, amongst oth- 
ers, in the high temperature superconductors, CMR ma- 
terials, quantum Hall bars, chemical and magnetic mix- 
tures, crumpled membranes, and some theories of struc- 
tural glasses [2] . Toy models in which these issues can be 
addressed are useful. In the current context, we choose 
these examples as they highlight subtleties typically ab- 
sent in the more standard spin models. In the continuum 
limit, any incommensurate rotationally symmetric model 
(with an interaction v(k) which is analytic about its min- 
ima and having non-vanishing second derivatives) will 
share much the same physics as the two specific models 
discussed below. When any incommensurate spin system 
possesses a rotational symmetry, the minimizing mani- 
fold of v(k) is a (d — 1) dimensional shell of radius q > 0. 
The Coulomb Frustrated Ferromagnet 
We now introduce the "Coulomb Frustrated Ferromag- 
net" . This is a toy model of a doped Mott insulator [4] , 
[6], of phase separation in High Landau level Quantum 
Hall systems [11], [12] and of certain amphiphilic sys- 
tems. It has been argued that within the Mott insulator 
the tendency, of holes, to phase separate at low doping 
is frustrated, in part, by electrostatic repulsion [4]. In 
three dimensions, a simple spin Hamiltonian [5] which 
represents these competing interactions is 

Hmou = - 2^ S(x)S(y) + — ^ 



(x,y) S^y 

= ^£[ A (*)+£i*-*r 2 ]|s(*)i 2 - 



(4) 



K 



In the second line we employed the Poisson summation 
formula. Here {K} is the set of all reciprocal lattice vec- 
tors as introduced in Section(II). Here, S(x) is a coarse 
grained scalar variable which represents the local density 
of mobile holes. Each site x represents a small region 
of space in which S(x) > 0, and S(x) < correspond 
to hole-rich and hole-poor phases respectively. In this 
Hamiltonian, the first "ferromagnetic" term represents 
the short-range (nearest-neighbor) tendency of the holes 
to phase-separate and form a hole-rich "metallic" phase, 
whereas the frustrating effect of the electrostatic repul- 
sion between holes is present in the second term. Non- 
linear terms in the full Hamiltonian typically fix the lo- 
cally preferred values of S(x). One may consider rf ^ 3 



dimensional variants wherein the spins lie on a hyper- 
cubic lattice, and the Coulomb kernel in Hq is replaced 
by — y\ 2 ~ d (or by ln|a; — y\ ] in two dimen- 

sions [15]) where = 2n d / 2 /T(d/2). Here the compe- 
tition between both terms, when Q < 1 favors states 
with wave-numbers ~ Q 1 / 4 . The introduction of the 
Coulomb interaction is manifestly non-perturbative: it is 
long range. Moreover, the previous ferromagnetic ground 
state becomes, tout a' coup, infinite in energy. We will, 
for the most part, focus on the continuum limit of this 
Hamiltonian where the kernel becomes 



' J Mott(k) 



Qfc- 2 + fc 2 (l+£^) 



(5) 



After rescaling, this may also be regarded as the small k 
limit of the more general 



v Q (k) = A(fc) + QlA(k)}- 1 + A[A(k)}< 



+AE(! 



coski)(l - cosfcj) + 0(k 6 ). 



(6) 



The constants A and A are pinned down if we identify 
v q{^) = VMott{k). Here, we will modify them in or- 
der to streamline the quintessential physics of this sys- 
tem. First, we set A = 0: In the continuum limit this 
term is not large nor docs it lift the "cubic rotational 
symmetry" [16] present to lower order. Next, we allow 
A to vary in order to turn on and off "cubic rotational 
symmetry" breaking effects [16]. The kernel vq{k) may 
be regarded as vm u augmented by all possible next to 
nearest neighbor interactions. As our Hamiltonian re- 
spects the hypercubic point symmetry group, by survey- 
ing all possible values of A wc should be able to make 
general statements regarding the possible phases (within 
the planes) of real doped Mott insulators. In dimension 
d = 2, whenever A > 0, the minimizing wave- vectors will 
lie along the cubic axis and "horizontal" order will be 
expected. When A < the minimizing wave-vectors lie 
along the principal diagonals and diagonal order is ex- 
pected. At large values of Q, when the continuum limit 
no longer applies, trivial extensions of these minimizing 
modes are encountered where one or more of the wave 
vector components is set to ir. The work reported in this 
paper focuses on spin systems on a lattice. The situa- 
tion becomes far richer in continuum field theories. In 
high level Quantum Hall systems, calculations on a sim- 
ilar model indicate [12] that the system undergoes a se- 
quence of transitions. At very low excess filling fractions, 
v — N < l/N (with N the integer part of the filling v) 1 a 
Wigner crystal of the cyclotron orbit centers is formed. 
As v is further increased, transitions between bubbles of 
increasing size occur. Ultimately, at a sufficiently large 
filling v, stripes appear. Similar textures are observed in 
chemical mixtures [2] which may be described by similar 
models [17]. We will prove, however, that this series of 
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FIG. 1. The shell of modes which minimize the energy in 
the continuum limit. For the case just discussed This sphere 
is of radius \q\ = Q 1 / 4 . 

transitions between various textures does not occur for 
multi-component spins on a lattice. On a lattice, barring 
special commensurate points that we detail below (e.g. 
the standard uniform and Neel phases and other, more 
intricate, commensurate states), only spiral (stripe-like) 
states may appear. 

If the ferromagnetic system were frustrated by a gen- 
eral long range kernel of the form V(\x — y\) ~ \x — y\~ p 
we could replace the [A(fc)] _1 ] in vq(k) by the more gen- 
eral [A(fc)]( p ~ d )/ 2 . For instance, in a spin-lattice ver- 
sion of the Quantum Hall problem, p = 1 and d = 2. 
Here, in the continuum limit, the minimizing modes are 
q ~ Q 1 /( 2 + d ~p) and as the reader will later be able eas- 
ily verify all our upcoming analysis can be reproduced 
for any generic long range frustrating interaction with 
identical conclusions. 

In Fig.(l), we schematically depict the manifold (M) 
of the minimizing modes in k space. When no symmetry 
breaking terms (A = 0) are present, in the continuum 
limit it is M is the surface of sphere of radius Q 1 / 4 . If 
A 7^ this degeneracy will be lifted: only a finite number 
of modes will minimize the energy. When A > there 
will be 2d minimizing modes (denoted by the big X in the 
figure) along the coordinate axes. In the up and coming 
we will focus mainly on A > 0. When A < 0, a mo- 
ment's reflection reveals that there will be 2 d minimizing 
modes along the diagonals, i.e. parallel to (±1,±1,±1) 
(and in this case, they will have a modulus which differs 
from Q 1/4 ). Unit ;ss explicitly stated otherwise, we will 
set A = for calculational convenience and when a finite 
A is invoked it will be made positive (to avoid the A de- 
pendence of \q\ incurred when the former is negative). At 
times, we will present results for vq(k) at sizable Q, even 
though the model was motivated as a good caricature of 
VMott only in the continuum limit (at small wave- vectors 

Membranes 

In several fluctuating membrane systems, the affinity 



of the molecular constituents (say A and B) for regions of 
different local curvature frustrates phase separation [2]. 
Let us define S(x) to be the difference between the A 
and B densities at x. In the continuum, the energy of 
the system contains a contribution, 

H mix = \j <?x \VS\ 2 (7) 

reflecting the demixing of A and B species. Instead of 
considering long-range interactions, we now allow for out- 
of-planc (bending) distortions of the sheet. Specifically, 
we assume that the two molecular constituents display 
an affinity for regions of different local curvature of the 
sheet. This tendency can be modeled by introducing a 
coupling term between the local composition S(x) and 
the curvature of the sheet. If the distortions remain 
small, we may write 

H c = J d 2 x [^a\Vh(x)\ 2 + |[V 2 /i(x)] 2 + AS(x)V 2 h{x)} 

= J d 2 xH c , 

where h(x) represents the height profile of the sheet (rel- 
ative to a flat reference state), a is its surface tension , 
and k is its bending modulus; A, the coefficient of the 
last term in the expression measures the strength of the 
coupling of the local curvature V 2 /i and the local compo- 
sition <fr, which we have included here to lowest (bilinear) 
order. This coupling term reflects the different affinities 
of the molecular constituents A (S = 1 corresponds to 
pure A composition) and B (S = corresponds to pure 
B composition) for, respectively, convex (V 2 h > 0) and 
concave (V 2 /i < 0) regions of the interface. The vari- 
ational Eqs. for the total energy H = + H c , with 
respect to the membrane shape {h(x)} read 

-crW 2 h + KW 2 (W 2 h)+AW 2 S = 0. (8) 

If |kV 2 (V 2 ^)| < min{|crV 2 /i|, lAV 2 ^} then 

AS{x) = ah(x) +g{x), (9) 

with g(x) a harmonic function satisfying boundary con- 
ditions. In Tl c , after an integration by parts, 

ASV 2 ^ w AS-V 2 S^ - — (VS) 2 . (10) 
a a 

H mix +H c k J d 2 x {h'\VS\ 2 + ^(V 2 S) 2 }, (11) 

where b' = b — —. This effective Hamiltonian, which is 
a function of S alone, reads 

H= I d 2 kv membrane (k)\S(k)\ 2 , (12) 
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where v mem brane(k) — \k 2 + |j#fc 4 is the 2D Fourier 
transform. A negative b' obtained when b < A 2 /a, sig- 
nals the onset of a curvature instability of the sheet. This 
instability generates a pattern of domains that differ in 
composition as well as in local curvature and thus as- 
sume convex or concave shapes. The characteristic do- 
main size corresponds to the existence of the minimum of 
the free energy at a non-zero wave number. The modu- 
lation length d — q~ x ~ y/(A 2 n/cr 2 )/\b'\ with q = \q\ the 
minimizing wavcnumber modulus of v mem brane(k) ■ After 
scaling, this model may be regarded as the continuum 
version of the frustrated short range kernel 

v z (k) = zA 2 (k)-A(k), (13) 

(where z — —A 2 n/{a 2 b')) on the lattice. The real-space 
lattice Laplacian 

All matrix elements {x\A R \y) = for \x — y\ > R = 
Range. Our system is of Range = 2. Explicitly, 

(x\A 2 \y) = 2d(2d + 2) for x = y 
— 4d for \x — y\ = 1 

2 for (x - y) = (±e t ± e e >) where I ^ (' 
1 for a ±2ee separation. (15) 

We will extend the investigation of this model over a 
broader range of parameters than suggested by its initial 
physical motivation. In the continuum limit, theories 
with high order derivative terms generally give rise to 
v(k) = P(k 2 ). where P is some polynomial. Although 
v z (k) and its likes are artificial on the lattice, their con- 
tinuum limit is quite generic. Later on we will show that 
if P(k 2 ) attains its global minima at finite |fc|, then ther- 
mal instabilities can incur extremely low values of T c . 

IV. ISING GROUND STATES 

In most sections to follow, the bulk of our reported 
results pertain to continuous spins. We now, however, 
present a short overview of discrete spin systems in or- 
der present new results and to illustrate how Ising sys- 
tems may be addressed within a similar momentum space 
framework. Towards the end of this section, will estab- 
lish the existence of a new "holographic" ground state 
entropy in certain incommensurate Ising systems wherein 
the ground state entropy scales as the surface area of the 
system. 

In an "Ising" system S(x) = ±1 at all lattice sites 
x. Stated alternatively, the scalar (n = 1) spins sat- 
isfy a normalization constraints {S 2 (x) — 1} at all N 
lattice sites x. Henceforth, we will adopt the latter 



point of view. The set of minimizing wave-vectors q, 
v(q € M) = mingli^fc)} defines a manifold M. If the 
local normalization constraints are swept aside then it is 
clear that the ground states are superpositions of sinu- 
soidal waves with wave- vectors q £ M. One would expect 
this to be true, in spirit, also in the highly constrained 
Ising case, if v(k) is sharply dipped at its global minima. 
If, in a non-rigorous setting, we "digitize" a particular 
plane wave S(x) = sign(cos(qi ■ x)) and compare it with 
the exact (numerical) ground state, then we will find en- 
couraging agreement in certain cases. For instance, this 
gives reasonable accord when H = H Mott [1,5,6]. This 
Hamiltonian (with some twists) was investigated in [4] on 
a square (d — 2) lattice. In the continuum limit (i.e. if 
the lattice is thrown away) we might naively anticipate a 
huge ground state degeneracy- a "digitized plane wave" 
for each wave vector q lying on the (d — 1) dimensional 
manifold {Mq : q 4 = Q} This large degeneracy might 
give rise to a loss of stability against thermal fluctuations. 
Striped phases (i.e. "digitized plane waves") were found 
in virtually all of the parameter range [1,5,6]. Only for 
a very small range of parameters were more complicated 
periodic structures found. 

An intuitive feeling can be gained by considering a one 
dimensional pattern such as 

+ + -- + + -- + + --... (16) 

This pattern is a pure mode S per i d=4 (x) — v / 2cos[ : |x — 
j]. With this observation at hand, a double checkerboard 
pattern such as 

++--++-- 
++--++-- 
--++--++ 

-- + + -- + + (17) 
extending in all directions in the plane is given by 
S(x) = 2co4- Xl - -} cos[-x 2 -^} = 
cos[|(xi + x 2 ) - |] + cos[|(xi - x 2 )]. (18) 

Such a 4x4x4 periodic pattern in three dimensions would 
include the eight modes i(±7r, ±7r. ± ir). This example 
illustrates an simple premise. A system composed of a 
periodic building block whose dimensions p\ x p 2 x ps, is 
given by 

S(x) = S , Pl (xi)5p 2 (x 2 )S'p3(x3). (19) 

If a configuration S p (x) contains the modes {fc™} 
with amplitudes {S p (k m )}, then Fourier transform- 
ing the periodic configuration S(x) one will find the 
modes (ifc™ 1 , ±fc™ 2 , ±fc™ 3 ) appearing with a weight ~ 
\S Pl (ki) x S P2 (k 2 ) x S P3 (k 3 )\ 2 . For high values of the pe- 
riods p, the weight gets scattered over a large set of wave- 
vectors. If v(k) has sharp minima, such states will not 
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be favored. The system will prefer to generate patterns 
s.t. in all directions i albeit one p, = 1 (or perhaps 2). 
For a pi x p 2 x P3 repetitive pattern, the discrete Fourier 
Transform will be nonzero for only Yii=i Pi values of k. 
This trivial observation suggests the phase diagram ob- 
tained by [6] in two dimensions. The intuition is obvious. 
We have derived [5] , rigorously the ground states in only 
several regions of its parameter space (those correspond- 
ing to ordering with half a reciprocal lattice vector), and 
on a few special surfaces (corresponding to ordering with 
a quarter of a reciprocal lattice vector). In all of these 
cases the Ising states may be expressed as superposi- 
tions of the lowest energy modes cxp[iq m ■ x\. Lately, 
a beautiful extension was carried out by [18]. Another 
nice work, focusing on scalar fields in the continuum is 
[19]. In the current publication, however, we are con- 
cerned only with spins on a lattice. We now ask whether 
commensurate lock-in is to be expected. The energy of 
the Ising "digitized plane wave" on an L x L x L lattice 
where q = (qi,0, 0) with q\ = 2ir/m, with even m, reads 



£ = ^£«(fc)|s(fc)| 5 



- E 

2 

m 



k 

v(k = ( 2 -^,0,0)) 
| exp[27rr//m] — 1| 2 



E 

j — 1,3,. ..,m— 1 



sin (7T7*/m) 



(20) 



The lowest energy state amongst all states of the form 
considered is a possible candidate for the ground state. 

Unlike the previous paragraphs and those to follow, 
we now try to present the reader with an intuitive feel- 
ing. The following paragraph is strictly non-rigorous. For 
the particular model long-range introduced in Eq.(6), it 
seems that for small values of Q, it might be worthwhile 
to have an incommensurate phase. This is, in a sense, 
obvious- all low energy modes are of very small wave- 
number q and hence not of low commensurability. In a 
stripe phase having of size 27r| q gr0 und the energy 



E 



-y OO 

N E 



n=0 



16v(k = (2n + l)q g round) 
(2n + 1) 2 tt 2 



(21) 



In the example of Eq.(6) for q gr0 und *C 1) the higher 
harmonics k = (2n + l)q, do not entail high energies. 
For large values of Q, q gr0 und — 0(1), and v(k = 
(2n + l)q gr0 und) can be very large if [(2n + l)q gr ound] 
approaches a reciprocal lattice vector K. Under these 
circumstances it will pay off to have a commensurate 
structure; for a u± x ui x ... x Ud repetitive block only the 
modes k = 27r(^, 2|, ^i) will be populated (i.e. have 

a non- vanishing IS^fc)! 2 )- the ferromagnetic point (a re- 
ciprocal lattice point) will not be approached arbitrarily 



close- if that is not true weight will be smeared over ener- 
getic modes. Generically, we will not be expect commen- 
surate lock-in in lim q — > for any theory with a frustrat- 
ing long range interaction. Although we have considered 
only striped phases (which have previously argued are 
the only ones generically expected), it is clear that this 
argument may be reproduced for more exotic configura- 
tions (such as bubbles, dots, cylinders etc.) As an aside, 
by replacing Fourier lattice sums by integrals, we note 
that in the continuum limit of Hmou, the ground state 
wave- vector in three dimensions, is immediately seen to 
be q gr0 und ~ Q 1 ^ 3 [20]- not the value suggested by dimen- 
sional analysis (if the lattice constant units are ignored) 
- the minimizing wavenumber q ~ Q 1 / 4 . We will prove 
that for continuous spins q gr0 und = 9- Similar results 
are obtained for the full Hamiltonian Hmou (riot only its 
continuum limit) and related systems, e.g. [18], [17], [21]. 

For finite range interactions it is easy to prove, by 
covering the system with large maximally overlapping 
blocks, that there will be a regime about q — 0, for which 
we will find the ferromagnetic ground state. A polyno- 
mial in A(fc) will have its minima at A(fc) = const, i.e. on 
a (d-1) dimensional hypersurface(s) in fc-space or at the 
(anti)fcrromagnetic point. The kernel v z (k) = zA 2 — A 
has its minima (z > 0) at q £ M z : A(q) = minj^j, 4d}. 
For z > gg: M z is (d — 1) dimensional. We may divide 
the lattice into all maximally over-lapping 5 x 5 x ... x 5 
hyper-cubes centered about each site of the lattice. 



Energy 



1 



5 x 6 d -! 



e(hypercube) (22) 



hypercubes 



and evaluate the energies e of all 5 x 5 x ... x 5 Ising 
configurations. Of all 2 5 configurations the Neel state 
will have the lowest energy for a sliver about z = 
Analogously for z > z top » 1, by explicit evaluation, 
the ground state will be ferromagnetic. Contour argu- 
ments can be employed and a finite lower bound on T c 
generated. 

If such a system satisfies periodic boundary conditions 
along the two square lattice diagonals e± : defined by 
xi T X2 = const, (with Xi the Cartesian coordinates), 
then an exponentially large number of Ising ground states 
can be constructed whenever z = |. To illustrate this, 
we note that when z = \, all minimizing modes lie on 



qe M z 



1 9i ± 92 1 



(23) 



By prescribing an arbitrary spin configuration along X- 
and fixing S(x-,x+) = S*(x_ , 0)(-l) x + : 



S(k) = S(x-,0) exp(ifc_x_) £(-lf + exp(i/c + .x + ) 

X- 35+ 

(24) 
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A) FLIP SPINS DIAGONALLY 




B) 

Q©0 

FIG. 2. Simple ground states for v z (k) with z = 1/8 (and 
for vq(Ic) when Q — 16): A) Their construction- along the ar- 
rowed line we arbitrarily prescribe spins. For each move along 
the dotted diagonal lines we flip the spins. B) In this chain 
of diagonal super-spins there is no stiffness against flipping. 

vanishes for |fc+| = \ki + k 2 \ 7^ tt- Similarly, by tak- 
ing the transpose of these configurations we can generate 
patterns having S(k) = unless |fc_| = \k\ — k 2 \ = n. 
The ground state degeneracy is bounded from below by, 
the number of independent spin configurations that can 
be fashioned along x + or x_, (2 i+1 — 2) where L is the 
length of the system along the x± axis. The number of q 
values, commensurate with the diagonal periodic bound- 
ary conditions, lying on M z= i is (4L — 2). We have thus 
proved that the ground state entropy in this two dimen- 
sional system is, at least, linear in the perimeter of the 
system. In a later section, we will later demonstrate that 
in the large n limit, such "holographic" like entropies do 
generically arise whenever the minimizing manifold M is 
(d-1) dimensional. The Ising ground state degeneracies 
are bounded from above by those of the spherical model 
(as any Ising configuration is a viable configuration in 
the spherical model). As the ground state entropy for 
our model is bounded both from above and from below 
by the perimeter of the system, the ground state entropy 
rigorously scales as the system size L. 

Similarly, in dimensions d > 2, we can set (d — 2) of 
the q components to zero. There are d(d — l)/2 cross- 
sections of the d-dimensional M,_\, all looking like the 

z — g o 

the two-dimensional M just discussed (i.e. |gi±g2| = 77 )• 
The real-space ground state degeneracy is bounded from 
below by d(d — 1)[2 L — 1] (along the (d — 2) zero-mode 
directions the ground state spin configurations display 
no flip). If we regard each diagonal row of spins as a 
"super-spin" then we will see that flipping any "super- 
spin" entails no energy cost. This is reminiscent to a 
nearest neighbor Ising chain where the energy cost for 
flipping a spin is dwarfed by comparison to the (logarith- 
mically) extensive entropy. We might expect that here, 
too, ordering might be somewhat inhibited. In two di- 
mensions, 

lim M z :q 2 = ^-, (25) 
the "average" number of allowed q e M z values -C 



O(L) (and similarly for the onset lim j_+ M z : (q — 

(±7r,±7r)) 2 = (g^ — j^))- For a hypercubic lattice of 
size Li x L2 x ... x Ld in d > 2 many discrete recipro- 
cal points will give rise to the same value of A(fc) <~ k 2 . 
The proof is trivial: if all Li = L, then the number of 
possible k 2 values is bounded by dL 2 , whereas there are 
L d k— values. Therefore, on "average" , the number of k 
points lying on M z , or more precisely lying lying the clos- 
est to M z , s.t. |A(fc) — Jj| is min, is, at least, 0(L d ~ 2 ). 

Of these, n-(ra-f) wave-vectors, with rn (and z) denoting 
the number of identical components (and the number of 
zero components) of a certain k\ nearest to M z , are re- 
lated to k\ by symmetry. As we have stated previously, 
in the continuum limit (q — ► 0) any short range kernel 
(including this one) will have a uniform (ferromagnetic) 
ground state. However the impossibility of constructing 
ground states that contain only "good" Fourier modes 
S(k € M) when M shrinks to a curved surface enclosing 
the origin is more general and will proved in the next sec- 
tion. For this short range model, even for z ^ g a huge 
ground state degeneracy is expected. A "plane wave" 
might correspond to each wave-vector q (or commensu- 
rate wave- vectors nearby) lying on the (d— 1) dimensional 
manifold M. As we shall prove later on, even in high di- 
mensions, and even if the interactions are long ranged, in 
the continuum limit it will not be possible to construct 
Ising states in which S(k £ M) = unless the minimizing 
manifold M contains flat non-curved segments (or more 
generally intersects a plane at many points). 

V. A UNIVERSAL PEIERLS BOUND 

In this section, we illustrate how a Peierls' bound can 
be obtained in a host of systems with two spin interac- 
tions. Although we specialize to translationally invariant 
two spin interaction in this section, as we will show ar- 
gue later, the technique that we outline below can be ex- 
tended to non-translationally invariant systems as well. 
In Peierls' elegant proof of long range order in the two and 
higher dimensional Ising models at sufficiently low tem- 
peratures, one of the necessary ingredients is the Peierls' 
bound. This bound amounts to the demonstration that 
having one ground state domain surrounded by another 
costs an energy which is, at least, linear in the length 
of their interface. In this section, we prove that such a 
bound is quite generic to many long range systems. Al- 
beit suggestive, the bound, on its own, does not suffice 
to demonstrate long range order. The unfamiliar reader 
may peruse [22]. If a real hcrmitian kernel v{k) attains its 
minima in only a finite number of commensurate recipro- 
cal lattice points {<&}, then a Peierls bound can, in some 
instances, be proven for an infinite range model: When 
possible this is suggestive of a finite T c . For instance, 
the bound for a (lattice) Coulomb gas (with the kernel 
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solving the discrete Laplace equation on the lattice) is 
trivially generated. 

v e (k) - v e (q) = e/A(k) - e/A(q = ir,ir,ir) 

>-A(A(k)-A(q = ir,n,ir))), (26) 

with A = 16d 2 . Here, the Hamiltonian 

H = ^E*(*W)i a - ^Exk'^ )|2 (27) 

depicts a lattice Coulomb interaction (wherein the spins 
portray charges). The right hand side of Eq.(26) is the 
kernel of an antiferromagnet. Both system share the 
same ground states. For a given configuration the en- 
ergy penalty for the Coulomb gas 

AE e = l/(2J\0$>e(fc) - Ve(q)}\S(k)\ 2 (28) 
k 

is bounded from below by the corresponding penalty in 
an antiferromagnet of strength A. In d = 2 the contour 
penalty of the antiferromagnet is 2^4|r|, (|T| = length of 
the contour T). A similar trick may frequently be em- 
ployed when the minimizing wave-vectors attain other 
commensurate values. It relies on comparison to a short 
range kernel for which a Peierls bound is trivial. This can 
be extended quite generally. All translationally invari- 
ant system with commensurate minimizing wave- vectors 
(q = (0,0,... 0), (ir, 0, . . . 0), . . .) at which the minimum 
of v(k) at k = q is quadratic may be bounded by a kernel 
of a system having nearest neighbor ferromagnetic and 
antifcrromagnetic bonds for which the Peierls bound is 
trivial (linear in the perimeter of the domain wall T). 

VI. ISING WEISS MEAN-FIELD THEORY 

We now estimate the critical temperature for all in- 
commensurate spin system in which within mean-field. 
If, in these systems, when T < T c , the on-site magneti- 
zation (S(x)) = s sign[cos(q ■ x)] (as suggested by some 
of the examples discussed hitherto) , then the sum 

J2(S(y))V(x = 0,y) = ^(S)(k)v(-k) 

V k 

°° 4(— I s !" 

= 3 E (2n+l)n [v{k = ^ + m) + <k = ^ + ml 

where we assumed q% = — with > 1 for all d com- 
ponents in replacing a discrete Fourier transform sum by 
an integral. The self-consistency equation 

s = (S(x = 0)> = -tanh[/?£V(* =0,y)(S(y))} (29) 



yields 

^ = \±^^v { k=(2n +m \. (30) 

n=0 v ' 

The non-trivial solution (with a non-zero magnetiza- 
tion s) should be self-consistent at all lattice sites (not 
only at x = 0). Self consistency at other values of x will 
lead to other lower bounds on T c which will read 

°° 8(— 1)™ 

1 E ( 9n , = {2n + W COS ^ 2n + A- (31) 



VII. 0(N > 2) GROUND STATES 

Henceforth, the bulk of the article focuses on continu- 
ous spin systems having n > 2 components. In an 0(n) 
spin system, the spins {S(x)} have n components and 
are all normalized to unity- S 2 (x) = 1. 

A common intuitive picture held by many is that in 
complex systems of the type presented in Section(III) is 
the complex ground states can arise. For instance, in 
real systems, nanoscale phase separation in the models 
of Section(III) can proceed through a certain class of vi- 
able ground states. Wigner crystals and arrays of dots 
are anticipated at large q. As q is decreased, the spatial 
size of the dots increases, gradually leading to a state 
approaching macroscopic phase separation. Phase dia- 
grams exhibiting devil staircases are not uncommon in 
scalar systems. Albeit their intuitive appeal and occur- 
rence in many continuum systems having scalar order pa- 
rameters [2], [5], [6], [23], ground states such as these do 
not generically occur in any 0(n > 2) system on a cubic 
lattice. In the up and coming we will prove that, barring 
special commensurability effects, any translationally in- 
variant 0(n > 2) system will generically display only 
spiral (stripe like) or screw like [24] ground states. For a 
discussion of the occurrence of dots, Wigner crystals etc. 
in the scalar case, the reader is referred to [2], [5], [6], [19]. 
Configurations such as these were detected numerically ( 
[4]) and analytically for the scalar Spherical model ( [5]) 
for the spin models that we discuss. We will explicitly 
prove, for the first time, that these do not generically oc- 
cur in incommensurate continuous finite n spin systems 
by deriving the relations to be satisfied by the ground 
state configurations containing a number of pairs of min- 
imizing wave-numbers {±<fi}. We will illustrate that spe- 
cial Wigner crystalline like patterns may occur for spe- 
cial commensurate wave-numbers alone. We explain how 
these can arise in the spherical limit in Section(XIV). 

We first note that our previous ansatz S(x) = 
sign(cos(qi ■ x)) is readily fortified in the 0(n > 2) sce- 
nario: here there is no need to "digitize"- in the spiral 
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Si(x) = cos(g • x), S 2 (x) = sin((f • x), S i>2 (x) = (32) 

the only non-zero Fourier components are S(q), S(—q). 
In plain terms, this state can be constructed with the 
minimizing wave- vectors only. It follows that any ground 
state g must be of the form 

S l (x) = J2< cos(g m -2+<t>T)- (33) 

m 

We now turn to the normalization of the spins at all sites, 

n 1 n 

l(x) = J2 S iW=2^^ aT <' 

i—1 m,m' i—1 

(cos(0™ + <f>f)cos[(q m +q m >)- A 
- s\n(4>? + 4>m') sin[(<f ro + q m >) ■ x] 
+ cos(0™ - (f>f)cos[(q m -q m >)- x] 
- sin(C - 4>f) sm[{q m - q m ,) ■ x\). (34) 

If l(x) = 1 is to hold identically for all sites x, then 
all non-zero Fourier components must vanish. For the 
{cos(A • x)} Fourier components: 

n 

°=[ E Y,a?a? coaW + t?) 

n 

+ E £ar<*'cos(C-C')L (35) 

and a similar relation is to be satisfied by the {sin(A-x)} 
components: 

n 

° = [ E E a >r'sm(c + c') 

n 

+ ^ ^ar< l 'sin(C-C')] (36) 

arcosC = <, arsinC^«r- (37) 
The {cos(A • x)} and {sin(A • a?)} conditions read 

n 

°= E E[<x"'-<x"'] 

n 

+ E £[«*'+ «r«r'] 

n 

°= E E^«r'+«r^] 

n 

+ E E^ m c'-c«r'] 

9m-9 m / = J i i=1 



Let us now consider the case of two pairs of minimiz- 
ing modes. For two pairs of wave-vectors ±<fi and 
±q 2 minimizing v(k), both not equal half a recipro- 
cal lattice vector: (0,0, ..,),(tt,0, ...,0), (0,tt,0, ...,0),..., 
(tt, it, 0, ...,0), (w, 7r, 7r), the vector A (up to an ir- 
relevant sign) may attain four non-zero values: A = 
2q\,2q 2 ,q\± q 2 . 

When A = qi + q 2 , the conditions are 

n n 
i=l i=l 

When A = <fi — (f 2 , these conditions read 

n n 
i=l i=l 

For,4 = 2<f Q (a = 1,2): 

n n 

= E t« - «] ,0 = 2^ « (40) 

i=l i=l 

Next, we define 



[/« = « =1 ,^,...,o 
i? a = K,..,<) 



(41) 



The previous conditions for A = <fi ±52, 2qi i2 imply that 



V 1 


■ u 2 


= u l 


■V 2 


= 


V 1 


■V 2 




■ u 2 


= 


u 1 - 


V 1 


= u 2 - 


V 2 - 


= 0. 



(42) 

The four vectors {U 1 , U 2 , V 1 , V 2 } are all mutually or- 
thogonal. The number of spin components n > 4. Two 
additional demands that follow are 

v a ■ v a — u a ■ Jj a 

2 

J2^° ■ V* + U a ■ U a ] = 2j2v a ■ V a = 2. (43) 

a—l a 

The last equation is the normalization condition- the 
statement that the coefficient of cos(A • x), when A = 0, 
is equal to 1. For the case of a single pair of wave- vectors, 
±<fi , A = 2qi , and the sole conditions are encapsulated 
in the last of equations(42) and in equation (43). A mo- 
ment's reflection reveals that this only allows for a spiral 
in the plane defined by U 1 and V 1 . When n < 4 there are 
no configurations which satisfy S 2 (x) = 1 identically for 
all sites x [excusing those having 2(qi+qj) = A is equal to 
a reciprocal lattice vector] that are a superposition of ex- 
actly two modes. For instance, with the convention that 
a "double checkerboard state" is a Neel state of 2 x 2 
blocks [6] (pi = p 2 = 2 in the two dimensional planar 
version of Eq.(19), we note that a (double checkerboard 
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state along the i = 1 axis) ® (a spiral in the 23 plane) 
has pairs (i, j) is a configuration in which A = 2(qi + qj) 
is a reciprocal lattice vector. As the number of minimiz- 
ing modes {q m } increases, some of the conditions may 
degenerate into one, e.g. if (qi + 5*2) = (93 — 9*2) (i-e. the 
modes are collinear). This degeneracy is the a second 
route that might allow for Ising configurations which are 
superpositions of several "good" minimum energy modes 
cxp(iq ■ x). The highly degenerate Ising ground states 
that we have constructed previously can fall under either 
one of these categories. If neither one of these situations 
occurs, Ising states cannot be superpositions of several 
minimum energy modes: we will be left with too many 
equations of constraints with too few degrees of freedom. 
For three pairs of minimizing modes, none of which is 
half a reciprocal lattice vector, {±q m }^ n=1 with 

q w ±q t ^q r ±q s ^ 2q w (44) 

for all w 7^ t, and r 7^ s, conditions similar to those 
that previously written for A = q\ ± ^, now hold for all 
(q w ± q t ). 

U a -Up = U 2 a 5 a , 

v a -v = vl s a , 

U a -V(} = (45) 

The relation U a ■ V a = (a = (3 in the last Eq above) is 
enforced by setting A = 2q a . Thus, when exactly three 
pairs of minimizing wave- vectors satisfying the equation 
are present, the vectors {U a , V a } define a 6-dimcnsional 
space, and hence n > 6. For p pairs of minimizing wave- 
vectors, n must be at least 2p— dimensional. This bound 
is saturated when S is (a spiral state in the 12 — plane) 
<g> (a spiral in the 34 — plane) (a spiral in the 2p — 

1, 2p plane), i.e. 

(ai cos(<fi • x + <f>i),ai sin(<f • x + (f>i), ... 
, a p cos(q p ■ x + <fi p ), a p sm(q p ■ x + <fi p )) (46) 

with Y^ P a =i a a = 1- When wave- vectors with q w ± q t — 
q r ± q s or q w ± q t = 2q r are present, pairs of conditions 
degenerate into single linear combinations. Thus far we 
assumed that for all i and j, A = 2(qi + qj) is not a recip- 
rocal lattice vector, s.t. sin(A • x) is not identically zero 
at all x G Z d . We term the such a p — 2 configuration 
a bi-spiral. It is simple to see by counting the number 
of degrees of freedom for n = 4, that the bi-spirals over- 
whelm states having only one mode ±<fi. This is a simple 
instance of a general trend: High p states are statistically 
preferred. Moreover, as we shall see later, they are more 
stable against thermal fluctuations. 

In summary, we outlined a way to determine all 
0(n > 2) ground states, whether commensurate or in- 
commensurate, for a given kernel V(x,y) — V(x — 
y). Whenever n > 2, any ground state configuration 



can be decomposed into Fourier components, S g (x) = 
Y^i^i {cosl'fi • x] + 6jsin[gi • x] j. The vectors ql are cho- 
sen from the set of wave vectors which minimize v(k). 
Here, the modulus \A4\ is the number of minimizing 
modes (the "measure" of the modes on the minimizing 
surface M. As long as these wave-vectors qi which mini- 
mize v(k) are "non-degenerate" , in the sense that the sum 
of any pair of wave vectors, ql ± <jj is not equal to the sum 
of any other pair of wave vectors, and "incommensurate" 
in the sense that for all i and j, 2(qi + qj) is not equal 
to a reciprocal lattice vector, then the that the condition 
[S 9 (x)] 2 = 1 can be satisfied only if \M\ < n/2. (In our 
toy model of the doped Mott insulator, these conditions 
are always satisfied for Q < 4.) Thus, for n < 3 only sim- 
ple spiral (|.M| = 1) ground-states are permitted, while 
for n — 4, a double spiral saturates the bound. Thus, 
generically, for 2 < n < 4 all ground states will be spi- 
rals containing only one mode. The reader should bear 
in mind that in the usual short range ferromagnetic case, 
the ground states are globally SO(n) symmetric and are 
labeled by only (n— 1) continuous parameters. Here, for 
each minimizing mode there are (2n— 3) continuous inter- 
nal degrees of freedom labeling all possible spiral ground 
states. For n > 2 this guarantees a much higher degen- 
eracy than that of the usual ferromagnetic ground state. 
If there are many minimizing modes (e.g. if the minimiz- 
ing manifold M were endowed with SO(d— 1) symmetry) 
then the ground state degeneracy is even larger! When 
n > 4, there are (generically) even many more ground 
states (poly-spirals). These poly-spiral states have a de- 
generacies larger than those of simple spiral. Their de- 
generacy 

g=p(2n-2p-l)\M\P, (47) 

where \M\ is the number of minimizing modes. We just 
proved that if frustrating interactions cause the ground 
states to be modulated then the associated ground state 
degeneracy (for n > 2) is much larger by comparison to 
the usual ferromagnetic ground states. 

VIII. SPIN STIFFNESS 

In this section, we focus for the sake of concreteness 
alone, on the Coulomb Frustrated Ferromagnet. A sim- 
ilar analysis, with different results, may be reproduced 
for all other interactions. When Q > in the in the con- 
tinuum (small k limit) of the Coulomb Frustrated Fer- 
romagnet of Eq.(6), with all 0(k 4 ) and higher terms ig- 
nored, the minimizing modes of VQ(k) lie on the surface 
of a sphere {Mq : q 2 = \/Q}. As Q — > 0, this surface 
Mq shrinks and shrinks yet is still a (d— 1) dimensional 
surface of a sphere. When Q = 0, the minimizing man- 
ifold evaporates into a single point q = 0. This sudden 
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change in the dimensionality has profound consequences. 
As we shall see shortly, it lends itself to suggest (quite 
strongly) that order is inhibited for a Heisenberg (n = 3) 
realization of our model. Before doing so, let us indeed 
convince ourselves, on an intuitive level, that the large 
degeneracy in k— space brought about by the frustration 
gives rise to a reduced spin stiffness. We now assume an 
ordered spiral state of momentum q on a cubic lattice 
(after all, we proved that in the general incommensurate 
case, these are the only ground states), and we examine 
the energy cost of varying the modulation wave- vector in 
a direction parallel and transverse to q. The exact ori- 
gin of the spiral ground is completely irrelevant (whether 
it is the the Coulomb Frustrated Model, the membrane 
model or other very different models)- all that matters is 
that the interaction kernel attains its minima at finite in- 
commensurate momenta \q\ which cover a spherical shell 
of radius q. For general rotationally symmetric incom- 
mensurate systems, Q may be regarded as a convenient 
shorthand, Q = q 4 . The upshot of the below, rather 
trivial, calculation is that rotationally invariant incom- 
mensurate continuous spin systems are very susceptible 
to transverse perturbations- a transverse twist may incur 
no energy penalty in the thermodynamic limit. 



A. Longitudinal Spin Stiffness 

We first examine the energy cost of a longitudinal twist 
to find that for highly incommensurate frustrated sys- 
tems with shells of minimizing modes, the energy cost 
for such a perturbation is just the same as in a nearest 
neighbor three dimensional XY ferromagnet. In a longi- 
tudinal twist 

S(x) = cos(—x + qxjei + sin(— x + qx)e 2 + SS 
Li hi 

= ei cos(fc ■ x) + e 2 sin(fc • x) + SS (48) 



B. Transverse Spin Stiffness 

When subjected to a transverse twist, the system re- 
sponds with only a quartic restoring potential. For a 
transverse twist, 

S(x) = cos[— — + qy]ei + sin[— — + qy]e 2 + SS (50) 

we find that as the difference of the energy kernels, v(k = 
qe 2 + ^fei) - v(q) = 7^7^, the energy gain, ignoring 



SS, is AE 



. This energy penalty for a transverse 



twist vanishes as L — ► 00 in d — 3- there is a complete loss 
of stiffness against transverse fluctuations. As N = L d , 
this energy penalty scales as O(y-). We note, in passing, 
that, to this order, this energy penalty matches with that 
of an XY chain, 



(51) 



An effective reduction in dimensionality seems to occur- 
the system effectively behaves as a one dimensional sys- 
tem when exposed to transverse perturbations. 

In the general case, for all incommensurate systems, 
this simple spin wave analysis yields a response to a twist 
8 with an effective kernel Ei ow (S) — v(q + 5) — v(q). For 
systems with a spherical shell of incommensurate mini- 
mizing modes, Ei ow ~ ^4±(5]^ for momentum deviations 
transverse to the minimizing q, while Ei ow ~ A\\5?, for 
deviations parallel to q. Such a fluctuation spectrum is 
reminiscent to that obtained in smectic liquid crystals. It 
will be noted, however, that at this stage of approxima- 
tion, Ei ow 7^ ^4±<5^ + j4||(5m, for general deviations <5: the 
precise fluctuation spectrum is much softer than that of a 
smectic liquid crystal (e.g. E\ ow vanishes for 5 connecting 
q with another vector on the minimizing manifold). 



with k = (^- + q)e\. The energy cost of this twisted state 
relative to the ground state is 

A£[{5(£)}] = ± - v(q)]\S(k>)\ 2 (49) 



Ignoring SS contributions, 

AE= Y Hk)- V (q ) ] = ^=[( T + qr 



-2i2 



8tt 2 A 



It is readily seen that this energy gain exactly coincides 
with the energy gain incurred by a uniform longitudinal 
twist in the three dimensional nearest neighbor ferromag- 
netic XY system, S(x) = cos[^]ei + sin[ 2 7 p]e 2 . 



IX. THERMAL FLUCTUATIONS OF ANY 
TRANSLATIONALLY INVARIANT XY MODEL 

In the previous section, we examined the response of 
an arbitrary incommensurate system having a shell of 
minimizing modes to various twists by simply examin- 
ing the energy difference between a ground state spiral 
and a spiral of another wave-number. We found that the 
system was very unstable to transverse twists. We now 
go one step further and completely treat the "soft-spin" 
version of the classical XY model with arbitrary incom- 
mensuration (not only those with a shell of minimizing 
modes). Towards the end of this section, we investigate 
what our general results imply in the case of rotationally 
symmetric incommensurate systems (which have a shell 
of minimizing modes) where we will find a link to smectic 
liquid crystal thermodynamics. 
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In this general theory, we include the non-linear inter- 
action Hi, 



Hsojt = H +u J2[S 2 (x) - l] 2 = H + Hi 



(52) 



with small u > 0, and forget about the normalization 
conditions |5(x)| = 1 at all lattice sites x. (The normal- 
ized, "hard-spin" , version can be viewed as the u — > oo 
limit of the soft-spin model.) As we proved previously, 
the only generic ground states (for both hard- and soft- 
spin models) when the spins S(x) have two (and also 
three) components are spirals 



^ground- 



state ^ 



cos(g-x); S 9 2 round - 



state (x) = sm(q- x). 



These are also the lowest energy eigenstates of H so f t . 
We will shortly expand H so f t about these ground states, 
keeping only the lowest order (quadratic) terms in the 
fluctuations SS. The quadratic term in {6Si(k)} stem- 
ming from H so ft is the bilinear j^(5S) + M(5S) where 

(SS)+ = (dSi(-ki),6S 2 (-ki),5Si(-k 2 )SS 2 (-k 2 ), 

SSi(-k 3 ), 6S 2 (-k 3 ), SSi(-k N ),6S 2 (-k N )) (53) 

and the matrix M reads 



/ 4 . 


. 1 


i 


. . • \ 


4 . 


i 


-1 . 




1 -i . 


. 4 


. 


1 i 


-i -1 . 


. 


4 . 


. i -1 




. 1 


—i 


. 4 


\'. '. '. 


. —i 


-1 . 


.0 4 / 



The sub-matrices are (2x2) matrices in the internal 
spin indices. The off diagonal blocks are separated 
from the diagonal ones by wave- vectors (±2q). Note 
that (fc|M|fc') = M(k - k'). Making a unitary (sym- 
metric Fourier) transformation to the real space basis: 
\x) = N^ 1 / 2 e lk ' x \k), the matrix M becomes block 
diagonal, (x\M\x') — M{x)5g^i. Diagonalizing in the 
internal spin basis, we find that 



6,2. 



(54) 



These eigenvalues may be regarded, in the usual ferro- 
magnetic case (the limit q = 0) as a two step (state) 
potential barrier separating the two polarizations. I.e., 
the normalization constraint of the XY spins (embod- 
ied in Hi) gives rise to an effective binding interaction. 
As we shall later see, when the number of spin compo- 
nents n is odd, one spin component will remain unpaired. 
Hi literally "couples" the spin polarizations. Employing 



Eq.(54), the corrected fluctuation spectrum {ipm}m=i (to 
quadratic order) satisfies a Dirac like equation 

[U+ v(-id x ) U + 2u ^ ] U+\i> m (x)) 

= E m U+\^ m {x)). (55) 



Here, 



U 



s'm(2q-x) 
2 cos(q-x) 
-1 — cos(2q-x) 
2 cos(q-x) 



1— cos(2g*-S") 
2 sin(q-x) 



(56) 



Alternatively, expanding in the fluctuations 5S(k), leads 
to bilinear (SS)+H(5S) where 



n k,k 



v(k) + 8u 






v(k) + 8u 



(57) 



along the diagonal, and Hj: ^ ±2 q = ^ u ( a 3 ^ ^i) on " 
the diagonal. Next, we perform a unitary transforma- 
tion exp(if CTi)Wg | g ±2 -.exp(-if cti) = 2u<r ± , while 
of Eq.(57) is unchanged. Till now, all that we stated, 
held for arbitrarily large u- our only error was neglecting 
0((6S) 3 ) terms by comparison to 0((SS) 2 ). Note that 
the main difficulty with the approach taken till now was 
the coupling between k and k ± 2q: i.e. k is coupled to 
k±2q, while k + 2qis coupled to k + Aq and k, and so on. 
Unless q is of low commensurability an exact solution to 
this problem is impossible. Equations with similar struc- 
ture appear in very different arenas (e.g. two dimensional 
Bloch electrons in a magnetic field [27], [28]) and such 
systems are often addressed via a continued fraction rep- 
resentation or by mapping the problem onto a Harper like 
Equation. To make rapid progress and to elucidate the 
similarity between incommensurate systems and smectic 
liquid crystals, here we will address our problem is a very 
direct and short fashion. Let us assume that u is small. 
In this case the lowest eigenstates of the fluctuation ma- 
trix will contain only a superposition of the low lying 
k- states (i.e. those close to the (d — 1) dimensional M 
(q > 0)). If ki = q + 5 is close to M, then the only 
important modes in the sequence {Si(k = ki + 2nq)} are 
ki, and k 2 = ki — 2q = —q + 5. The sub-matrix in the 
relevant sector reads 

/v(ki)+8u Au \ 

v{ki) + 8u 

v(k 2 ) + 8u 

V Au v{k 2 ) + 8uj 

The lowest eigenvalue reads 



E iow = -Hki) +v(k 2 )] +Au 
-l\/[v(ki)-v(k 2 )} 2 +6Au 2 . 



(58) 
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Equivalently, this can be determined from the direct com- 
putation of the determinant to 0(u 2 ): to obtain 0(u 2 ) 
contributions we need to swerve off the diagonal twice. 

N 

dctH = Y[[v(ki) + 8u} 2 - (4w) 2 ^ [v(kj) + 8u] 

i=l j 

x [v(kj + 2q) + 8u] [j [v&i) + 8u} 2 

ki^kj ,kj+2q 

The fluctuation spectrum is trivially determined by re- 
placing v(k) by [v(k) — E] in det Ti. and setting it to zero. 
To this order we re-derive Ei ow . Higher order terms in 
the determinant may be trivially computed. The parti- 
tion function is Z = const [detH]^ 1 ^ 2 . For any u, no 
matter how small, there exists a neighborhood of wave- 
vectors k near q such that \v(k) — v(k + 2q) \ <C u and as 
before we may re-expand the characteristic equation for 
these low lying modes, solve a simple quadratic equation, 
expand in the components of (k — q) and obtain a simple 
dispersion relation. 

Till now, our results, held for arbitrary incommensu- 
rate XY systems. We now investigate what transpires 
when the minimizing modes form a continuous shell M 
(e.g. all incommensurate systems with a rotational sym- 
metry that secures, as q ^ 0, a shell of radius q of min- 
imizing modes about the origin). For concreteness, if 
5 = S± + 5\\ en, with _L , || denoting directions orthog- 
onal, and parallel tonl M, then by expanding Eq.(58) 
to lowest orders for the example presented in Eq.(6) 

E low = A ± 6i + A ll Sf ] , (59) 

where 

l d 2 v{\k\) A|| 

A " - 2^fc^ l|fe| ^' A± = W (60) 

Examining the fluctuations about an ordered state, we 
find that the effect of thermal fluctuations is reduced 
by comparison to the dispersion relation in the large n 
model and spin wave stiffness analysis. Nonetheless, at 
this level of analysis, order is still inhibited, much un- 
like many models wherein "order out of disorder" oc- 
curs [25] by entropic stabilization about certain viable 
ground states. This dispersion is akin to the fluctuation 
spectrum of the smectic liquid crystals [44] which is well 
known to give rise to algebraic decay of correlations at 
low temperatures. In our case, q ^ and thus the corre- 
lations should have an oscillatory prefactor. To be more 
concrete we find that, for general rotationally symmetric 
incommensurate spin systems having a minimizing shell 
of modes at |g| — q > 0, the correlator (in cylindrical 
coordinates) is 

Arp' x 2 O 

G{x) ee (S(0) ■ S(x)) ~ -5- expI-2^ - ^i(^)] 

x cos[gx||], (61) 



where r\ — ^§-q, d = ^ where A is the ultra violet 
momentum cutoff, 7 is Euler's const., 

00 „ 

El (z) = - 1 -lnz-J2(-r^- r (62) 

n=l V '* 

is the exponential integral. For the uninitiated reader, 
we present this standard derivation in the appendix. 

The thermal fluctuations / d d k/ Ei ow {k) diverge as 
[— In |e|] with a e a lower cutoff on \k — <j| [29] . Such a log- 
arithmic divergence is also encountered in ferromagnetic 
two dimensional 0{2) system if it were exposed to the 
same analysis. Indeed, both models share similar char- 
acteristics albeit having different physical dimensionality. 
Thus far, our analysis was performed on incommensurate 
structures with minimizing shells in fc-space. For frus- 
trated systems having special commensurabilities, intri- 
cacies may arise (some are explicitly detailed in [26]). 
Relations to Smectic Liquid Crystals 
We now fortify our earlier conclusion regarding a tight 
relation between incommensurate continuous spin sys- 
tems and smectic liquid crystals. For incommensurate 
spin systems with minima on a (d— 1) dimensional shell 
in k space, the resulting spectra attained by both the 
Dirac like equation and the simple more naive spin stiff- 
ness analysis, matches that of smectic liquid crystals. 
Smectic ordering involves a breaking of rotational sym- 
metry and of translational symmetry along one direction 
alone. Let u be the displacement of the smectic along 
the z-direction. When the wave-vector k is along the z- 
direction the displacement is longitudinal, and the energy 
is of the elastic form ^Bk 2 \u(k)\ 2 , where B is the com- 
pressibility for the smectic layers. When k is normal to 
z, the displacement is transverse to the layer separation. 
No second order in k± the displacement costs no energy. 
Here, the restoring force is associated with a director 
splay distortion leading to an clastic energy -f-(^ ' "-) 2 
with K the splay constant. In, e.g. Smectic A, as h is 
normal to the layers, Sn = — V^u and an elastic energy 
^Kk±\u(k)\ 2 . Thus in smectic liquid crystals, the kernel 

v(k) = Bkf\ + Kki. (63) 

This is much alike our dispersion for incommensurate 
phases apart from a shift: in the incommensurate phases 
the deviations were about a finite wavenumber q whereas 
in liquid crystal, the uniform k — is the ground state. 
As seen in the form that we obtained for G(x), the pen- 
etration depth A = \jKjB determines the decay of an 
undulation distortion (splay director distortion) imposed 
at the surface of the smectic. 

In conclusion, if H so f t is indeed soft (u <C 1) then, 
in rotationally symmetric incommensurate soft spin XY 
systems, quasi long range (algebraic) order may be ob- 
served at low temperatures, These systems essentially 
may smectic thermodynamics. 
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X. A GENERALIZED 
MERMIN-WAGNER-COLEMAN THEOREM 

We will now generalize the Mermin-Wagner-Coleman 
theorem [30], [31]: All continuous spin systems with 
translationally invariant two-spin interactions in two di- 
mensions with a rotationally symmetric twice differen- 
tiable [33] Fourier transformed kernel v(k) show no spon- 
taneous symmetry breaking at finite temperatures T > 0. 
This is true for all systems, irrespective of the range of 
the interaction or of its nature. As will be made clear 
shortly, our analysis holds for both commensurate and 
incommensurate systems. 



A. The Classical Case 

Our approach is the standard one. We will keep it more 
general instead of specializing to anti/fcrromagnetic or- 
der or to interactions of one special sort. With the nota- 
tions introduced in Section(II), we investigate n compo- 
nent spins on a lattice. An applied magnetic field 



h(x) = hcos(q ■ x)e c 



(64) 



causes the spins to take on their ground state values. 

If n = a = 2 the unique spiral ground state (S 9 (x)), to 
which a low temperature system would collapse to under 
the influence of such a perturbation is 

Sf(x) =sin(q-x). (65) 
When n = 3 the ground state is not unique: 
S !<n(x) = nsin(q-x) 

Y, r 2 i = 1 (66) 



and a magnetic field may be applied along two directions, 
with all the ensuing steps trivially modified. With the 
magnetic field applied 

H = \ E E v{n-mwsM - E MWf). (67) 

x,y i—1 x 

Note that the knowledge of the ground state is not im- 
perative in providing the forthcoming proof [32] (whether 
a spiral or any other state). Our analysis thus holds for 
both commensurate and incommensurate systems. We 
exploit the standard rotational invariance of the measure 

Jdfi ■= Z- 1 J ]Jd n S(x)5(S 2 (x) - l)e~ 0H . (68) 



This is not applicable to many other spin models (e.g. 
Dzyaloshinskii-Moriya interactions, orbital Jahn- Teller 



and Kugel-Khoskii models). The generators of rotation 
in the [a/3] plane at a lattice site x are 



Lf = S a {g) 



d 



dSP(x) 



S {x) 



d 



dS a (x) 



(69) 



For any single spin, 

f(Si, S a cos + Sp sin 6, ... 
, Sp cos 6 — S a sin a, S n ). 

It follows that 



(70) 



= J d n SS(S 2 -l)L af3 f(S). (71) 

In the up and coming, _L will denote the the projection 
along the (3 direction. We define the operators 

A(k) = ^ exp[ifc ■ x\S±(x), 

X 

B(k) = ^exp[i(fc + q) ■ x]L s ((3H), (72) 

X 

where L x = {Li}^ 1 with Li = L%= n ^=\ 
By the Schwarz inequality, 

K E A * B ^ 2 ^ ( E A i A ^ * ( E B*Bi). (73) 

i— a,/3 i=a,j3 i=a,(3 

We will let i = a, (3 in the sum span a two element subset 
of the n spin components. For any functional C: 

L s (e- fm C) = e-" H {L s (C) + CL S (-(3H)}. (74) 

Invoking Eq.(71), 

= /il dn ' S ^ 6 (S 2 (x) - 1) L s [e- fiH C\. (75) 



(76) 



Employing the last two equations in tow, 
(CB(p)) = (^exp[^-f] L S (C)). 

X 

It is readily seen that 

£<4(4(/?i0)) =/?([E Si^SMVix-y) 

i—a,(3 i—a,(3 

-h(x)S n (x)}), (77) 

and 

(B(k)* ■ B(k)) = /3^{(cos(fc + q) ■ (x - y) - 1) 

x,y 

[ Y (Si{mm]v(z-m 



i— a,/3 



-h{x)(S n {x)) > (78) 
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Henceforth, for simplicity, we specialize to n = 2. 
Fourier expanding the interaction kernel 



(79) 



and substituting 

(S(x) ■ S(y)) = -L^(\S(u)\ 2 )e iS < s -V, (80) 



we obtain 



< (B(k)* ■ B(k)) = pAfE - phfiSni-q)) = 



— [v(u + k)+v(u - k) 

U 

2v(u)](\S(u)\ 2 } - ^(S n (-q)) (81) 



(2) 

where A - E measures the finite difference of the internal 

k 

energy with respect to a boost of momentum k. 

x exp[zfc- (x - y)}. (82) 

(A(k)* ■ B{k)) = LUS±(S)) exp[i(fc + q) ■ £]) = m ?) 

where = (S n (q)} and, as noted earlier, _L refers to 
the i — 1 spin direction orthogonal to i = n = 2. Note 
that with our convention for the Fourier transformations, 
a macroscopically modulated state of wave-vector q, the 
magnetization m q = 0(N) as is the energy difference in 
Eq.(81). Trivially rewriting the Schwarz inequality and 
summing over all momenta fc, 



f(B(ky-B(k)) ~ Y 

which explicitly reads 

2N\m q f((3j2((\S(u)\ 2 )Hk + u) 

fc 

+v(u - k) - 2v(u)] + 2\h\\m q \ 



(83) 



k x >v 



Nj2(Sl(x)). 



(84) 



Explicitly, as the integral j' fe ' ><5 ... is non- negative 
(as (B(k)* ■ B(k)) > the denominator in Eq.(84) is 



positive for each individual value of k), and as (Sj_(x)} < 
1, we obtain in the thermodynamic limit 



\k\<6 d d k 



d d i 



[2-kY V J {2ttY 

(\S{u)\ 2 )(v(k + u)+ v(k -it)- 2v(u)) 

-l 

+2|ft||m,| < 1. 



(85) 



Taking S to be small we may bound from above (for each 
value of fc) the positive denominator in the square brack- 
ets and consequently 



/ 



\k\<5 d d k 



(2tt) c 



r d d u 

J wr 



(\S(u)\ 2 )(v(k + u)+ v(k -u)- 2v(u)) 



+2\h\\i 



> 



\k\<s d d k 



(2ir) d 



A 1 k 2 \ a {\S{u)\ 2 ) + 2\h\\m q 



d d u 
{2ir) d 
l 



(86) 



with Aa chosen to be the largest principal eigenvalue of 
the d x d matrix didj[v(it)], and Ai a constant. For a 
twice differentiable v(u), and for |fc| < 5 where 5 is finite, 

(v(k + u)+ v(k - u) - 2v(uj) < A x \ z k 2 < B x k 2 (87) 

for all u within the Brillouin Zone with the additional 
positive constant B\ introduced [33]. Here we reiterate 
that (B*(k) ■ B{k)) of Eq.(81) is positive definite and 
consequently the bound derived is powerful. 

In d < 2, the integral j' fe ' <<5 jL^ diverges making it 
possible to satisfy Eq.(85), at finite temperatures, when 
the external magnetic field h — > only if the magne- 
tization m q = 0. If finite size effects are restored, in 
a system of size N = L x L where the infrared cutoff 
in the integral is 0{ 2 ^, 2 f-) the latter integral diverges 
as 0(hiN). This implies that the upper bound on \m^\ 
scales as 0(N/^/\nN), much lower than the O(N) req- 
uisite for finite on-site magnetization. For further details 
see [32]. If there are M > 2 pairs of minimizing modes 
and 2p + 1 > n > 2p (with an integer p) then we may 
apply an infinitesimal symmetry breaking magnetic field 
along, at most, min{p, M} independent spin directions 
(a). Employing the spin rotational invariance within 
each plane [a(3] associated with any individual mode I 
(au lieu of a specific q) we may produce a bound simi- 
lar that in Eq.(85) wherein (|5'('u)| 2 ) will be replaced by 
E^d^)! 2 ) and |m 9 -| -> |m Q (/)|. 



B. The Quantum Case 

The finite temperature behavior of a quantum system 
is, in many respects, similar to that of a classical system. 
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The quantum system is also invariant under rotations 
with 5 2 (x) = 5(5 + 1). Alternatively, one could directly 
tackle the n = 3 quantum case by applying the Bogoli- 
ubov inequality 



({A,Al})*([ [C,^],Ct]> > |<[C,^]>| S 



(88) 



with [ , ] and { , } the commutator and anticommu- 
tator respectively. From this inequality, it follows that 
([ [C, H], C^]} is positive definite. In particular, for any 
six operators {Ai,A 2 , A 3 , C\, C 2 , C 3 }, 

§E<{4..4»)(E<[ [Ca,H],Cl\))> 



o=l 



^|([a,A 4 ])| 2 (89) 

a 

Setting A\ — S 2 (q — k) and C\ — S\(k) we will once 
again obtain Eq(85) with the classical spins replaced by 
their quantum counterparts. 

Rather explicitly employing 

[S a {k),SP(k')\ = ie a ^S 7 (k + k'), (90) 
we find for the Hamiltonian of Eq.(67) (with n = 3), 

[[CuH],C\] = ±^( S ^') S 2(-k') + Ss(k')S 3 (-k')) 



x[v{k + k') - 2v{k) + v{k' - k)] 

+ ^E /i 3(fc')^3(-fc')- 



Similarly 



Kld,^])!^!^^)! 2 , 



(91) 



the squared magnetization along the z (or 3) direction 
for a mode q, and 



^{{A U A\}) =2^(S 2 (k)S 2 (-k)) 
k k 

= 2([S 2 (x = 0)} 2 ). 



(92) 



Next, let us cyclically set, A 2 = Ss(q — k),C 2 = 
S 2 (k),A 3 = Si(q-k), and C 3 = S 3 (k). The com- 
mutators associated with these operators are all iden- 
tically the same apart from a uniform cyclic permuta- 
tion of all spin components involved. Trivially rewriting 
the symmetrized Bogoliubov inequality Eq. (89) and sum- 
ming over all modes fc, 



4 ti\ > y Eg l([Cq^a])|- 



k.a 



?Z a ([[Ca,H],Cl]) 



(93) 



Replacing the k sums by integrals in the thermodynamic 
limit, and employing the positivity of ([[C a , H], C^]) that 
follows from the Bogoliubov inequality for each individ- 
ual value of k, we find the quantum analogue of Eq.(85), 



*\k\<6 d d k 



2 /- |K|<0 d d k r f d d u 
2/3'"^ J J^YJ WF 
(\S{u)\ 2 )(v(k + u)+ v(k - u) - 2v{u)) 



+ 



m„ 



< 5(5 + 1). 



(94) 



Apart from the simple scaling factor of 45(5 + 1) = 
2 A^}} by comparison to the classical case, there 

is no difference between this inequality and its classi- 
cal counterpart in the zero field limit. We symmetrized 
the Bogoliubov inequality (Eq.(89)) in order to avoid the 
appearance of only two transverse spin components in 
the Fourier weights {\Si(q)\ 2 ) so as to give the resulting 
Mermin- Wagner inequality a transparent physical mean- 

(2) 

ing associated with the energy boost differences AL 'E 
which are symmetric in all spin indices. From here on- 
ward the discussion can proceed as in the classical case. 



XI. MERMIN- WAGNER BOUNDS IN HIGH 
DIMENSIONS 

In any dimension, Eq.(85) reads in the limit h 
d d k 



2m, z T 



f d d k 1 

J (2^ ^) 



< 1, 



E 



(95) 



with the shorthand defined by Eq.(81). By parity invari- 
ance and noting that {S(x)} are real, 

E v(k + u)(\S(u)\ 2 ) = E «(* " "Xl^)! 2 )- (96) 



Thus the denominator of Eq.(95) reads A^E = [2{E % - 
Eo)] where E^ is the internal energy of system after un- 
dergoing a boost of momentum k and Eo denotes the 
internal energy of the un-boosted system. In some in- 
stances when the dispersion relation about an assumed 
zero temperature ground state is inserted into Eq.(95), 
we will find that the integral in Eq.(95) diverges: At 
arbitrarily low temperatures we cannot assume the zero 
temperature ground state with the natural dispersion re- 
lation AE for fluctuations about it. A case in point is the 
dispersion relation for the Coulomb Frustrated Ferromag- 
nct. The denominator in Eq.(85) is a finite temperature 
extension of the T = dispersion relations (Al ' E). In 
general, at zero temperature, 



<|sW> 



N 2 



(97) 



1G 



and the integral in Eq.(95) becomes 



/ 



d d k 



1 



(27r) d v(k + q) + v(k - q) - 2v(q) ' 



(98) 



For the incommensurate models, the reader will recog- 
nize this as none other than the integral obtained by 
other methods via the Dirac like analysis in Section(IX). 
Specifically, the denominator of Eq.(98) is E [ow of 
Eq.(58) wherein ki = q + k and ki = —q + k (which 
transformed into Eq.(59) (with S now portrayed by k) 
for general rotationally symmetric incommensurate sys- 
tems). 

Whenever Eq.(98) diverges, an assumption of an al- 
most ordered ground state at arbitrarily low tempera- 
tures (T = + ) with the ensuing zero temperature dis- 

(2) 

persion relation (AVE = [2(E% — E )]) about that or- 
dered state is flawed: Eq.(95) is strongly violated. This 
poses no problem for most canonical d > 2 dimensional 
models where the integral in Eq.(98) is finite. For all 
three dimensional incommensurate systems with a rota- 
tional symmetry that secures, as q ^ 0, a shell of ra- 
dius q of minimizing modes about the origin (e.g. the 
three dimensional Coulomb Frustrated Ferromagnet and 
other high dimensional systems), the divergence of this 
integral hints possible non-trivialities. The decoherencc 
time scale (or more precisely, bandwidth time scale [34]) 
(t) may diverge if we assume a zero temperature disper- 
sion of fluctuations about an ideal crystal. A divergent 
decoherence time scale suggests the absence of broken 
translational symmetry. Formally, 



/ 



d 3 k 



r(r = 0). 



(99) 



In order for magnetization m q = O(N) to arise, the av- 
erage of the inverse boost energy over all of k space 



A ( - 2) E 



(100) 



at all < T < T c . In most d > 2 systems this is trivially 
satisfied with the average bounded by a constant at zero 
temperature. This average diverges at T = whenever 
the integral of Eq.(98) does. There are two possibilities: 

(i) The system is ordered at all temperatures T < T c 
in which case, the thermodynamic average of Eq.(lOO) is 
finite for all T > and r is non- analytic at T = 0. 

(ii) The system is disordered at all finite temperatures 
and orders classically only at T = 0. 

The first possibility ((i)) was argued for by a non- 
rigorous yet elegant diagrammatic analysis by Brazovskii 
[36] long ago: Thermal fluctuations, on their own, may 
enhance (or generate) cubic terms fortifying (or trigger- 
ing weak) first order transitions. This cannot be ruled 
out by the rigorous Mermin- Wagner inequalities that we 



derived. As reiterated, if the fluctuation integral diverges 
then we may not obtain the low temperature T = + 
dispersion by assuming a nearly perfectly ordered state. 
This does not rigorously preclude order. Order, if it exists 
at arbitrarily low temperatures must display non-trivial 
excitation spectra about it having a vital explicit thermal 
dependence. If such a possibility arises, it might be im- 
material if the system is permanently frozen into a glass 
before reaching an equilibrium thermodynamic transition 
[13]. In both cases, r(T) is non-analytic at T = 0. The 
integral of Eqs.(98,100) has a suggestive physical inter- 
pretation. If the quantum spin system is subjected to a 
boost of momentum k, then 



(101) 



is the characteristic lifetime of the excited state. The 
average in Eqs.(100) is the characteristic relaxation (or 
decoherence) time of the system averaged over magnons 
of all possible momenta. Whenever the average charac- 
teristic relaxation time 



f d d k 



d d u 
{2ir) d 



{\S(u)\ 2 }(v(k + u) + v(k -u)- 2v{u)) 
diverges then by our generalized inequality 



2\m q \ z T{Tr) = 2\m q \ z TT r=0 < 1 



(102) 



(103) 



the system does not order in such a way that the fluctu- 
ation dispersion about any viable ground state is valid. 
We emphasize that, even in case (ii), the Mermin- Wagner 
inequality allows for algebraic long range order. Indeed, 
we have found by our primitive n = 2 analysis (assuming 
a fluctuation spectrum about a T — ordered state), that 
although the average magnetization vanishes at all T > 0, 
the correlations decay algebraically (see Eq.(61) and the 
appendix where this correlator is explicitly evaluated). 
As fortold, the main difference between the results of our 
rigorous analysis here and that of the soft spin model 
in section(IX) is that the soft spin analysis relied on a 
comparison to the ground state so that no explicit ther- 
mal dependence of the mode occupancies (IS'(fc)l 2 ) was 
allowed. As fortold, if the mode occupancies have a non- 
trivial thermal dependence which cannot be neglected 
even at arbitrarily low T = + temperatures, ordering 
may occur. 



XII. CLASSICAL 0(N = 3) FLUCTUATIONS 

We now return to the more naive "soft-spin" (see sec- 
tion (IX) for the soft spin XY) 0(3) models in order 
to witness an intriguing even-odd binding-unbinding (or 
pairing - unpairing) effect that could have otherwise been 
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missed. We will see that Hi plays the role of a pair- 
ing interaction. We will see the incommensurate clas- 
sical even n soft-spin systems are more "gapped" than 
their odd n counterparts. When, in three dimensions, 
incommensurate continuous spin systems are endowed 
with a spatial rotational symmetry which secures a finite 
shell of minimizing modes, odd n systems may be dis- 
ordered while their even n counterparts may exhibit, by 
the same soft-spin analysis, critical algebraic long range 
order. This even/odd alteration is reminiscent of the 
appearance/non-appearance of a gap for quantum spin s 
chains with even/odd values of (2s + 1) [37], [38] and the 
physics for an even/odd number of legs in an s = 1/2 
quantum spin ladder [39] where such even/odd discrep- 
ancies are triggered by Berry phase terms. As before, 
we will begin our analysis for general incommensurate 
soft spin models and only at the end investigate what 
occurs for rotationally symmetric systems with a shell 
of minimizing modes. As we proved, for an n — 3 sys- 
tem the generic ground states are simple spirals. If we 
rotate the helical ground-state to the 1 — 2 plane; the 
single quadratic term in 5Si=3(x) is ^2s(SSi—3(x)) 2 . The 
eigenvalues Ai> 3 = 2 = A_ = A m j n (see Eq.(54). Here, all 
that follows holds for arbitrarily large u- the only approx- 
imation that we are making is neglecting 0((5S) 3 ) terms 
by comparison to quadratic terms: i.e. assuming that 
SS(x) <C 1. Unlike the above treatment of the XY spins, 
no small u is necessary in order to make headway on the 
Heisenberg problem. For n = 3 we find that the fluctua- 
tion eigenstates of the are the products of an eigenstate of 
H within the plane of the spiral ground state and a fluc- 
tuation eigenstate of v(k) along the direction orthogonal 
to the spiral plane. Written formally, to quadratic order, 
the fluctuation eigenstates are \ip m ) <8> \SSs(k)). Fluctu- 
ations along the i = 3 axis are orthogonal (in a geomet- 
rical and formal sense) to the ground-state plane. This 
is expected as fluctuations in any hyper-plane perpen- 
dicular to the [12] plane do not change, to lowest order, 
the norm of the spin. |<55 3 (fc)) is literally the "odd" man 
out. As foretold, this is a general occurrence. Whenever 
the number of spin components is odd, one unpaired spin 
component is unaffected by the interaction enforcing the 
spin normalization constraint. 

Within the i = 3 subspace (k\M\k') — 2Sj: p and our 

previous analysis follows. The dispersion EfT 3 = [v(k) + 
2 u Xi—s] does not have a higher minimum than with E m 
[40]. Both Ejr 3 and E m share the same value of A and 
consequently the SSi = s(x) fluctuations are minimized at 
wave- vectors I G Mq such that v(£) = v(q) = mm^v(k). 
As IVVn) is a normalized superposition of \SS(k)} modes 
(the latter spin vectors being in the 1 — 2 plane) and v(k) 
is diagonal, and attains its minimum at £ € Mq: 

rmn{E m }>E^ MQ . (104) 
Here, we invoked the trivial inequality 



min^lItfo + tfiM >mm<0|JT o |0>+niin<£|iri|O- (105) 

4> <t> ? 

Thus, there exist Goldstone modes corresponding to 
SSi=3 fluctuations, and one must adjust additive con- 
stants such that minr{iv =3 } = 0. The resulting fluctua- 
tion integral reads 

/d 3 k 1 
T^W "7a 7v (106) 

where we invoke (A^(fc)A^(fc')) = S U p fi (\ASi(k)\ 2 ) 
by translational invariance, and employed equipartition. 
Till now we investigated general incommensurate sys- 
tems, we now specialize to rotationally symmetric sys- 
tems. As pointed out earlier, when the incommensurate 
spin system possesses a rotational symmetry, the mini- 
mizing manifold M is a (d — 1) dimensional shell of ra- 
dius q > 0. The fluctuation integral receives divergent 
contributions from the low energy modes nearby. By 
quadratic expansion about the minimum along n _L M, 
a divergent one-dimensional integral for the bounded 
((AS(x = 0)) 2 ) signals that are quadratic fluctuation 
analysis calculation is inconsistent. We are led to the 
conclusion that higher order constraining terms are im- 
perative: We cannot throw away cubic and quartic spin 
fluctuation terms (SS 3 ' 4 (x)) relative to the quadratic 
(5S 2 (x)) terms (notwithstanding the fact that all of these 
terms appear with 0(u) prefactors irrespective of how 
large u is). The spin fluctuations 5Si = z{x) are of order 
unity at all finite temperatures and T c (q > 0) ~ 0. Note 
that this is "almost a theorem" . Here we do not demand 
that u be small (only 5S). This is an important point. 
\ip m ) is an eigenstate for arbitrarily large u. To quadratic 
order in the fluctuations, the minimum belongs to \8S?,) 
(or is degenerate with it). The divergent fluctuations 
here signal that 0(SS 4 ) = 0(5S 2 ). Thus assuming that 
SS < \fjju (with J = 1 the exchange constant) we reach 
a paradox. Thus, if the integral in Eq.(106) d iverg es for 
an 0(3) model then at all finite T: SS > \J J/u. Un- 
like the case for n = 2, there is no algebraic long range 
order. Computing the correlator, we find that G(x) de- 
cays exponentially at all temperatures (much unlike the 
algebraic long range order seen in Eq.(61)). 

When the incommensurability q vanishes, the minimiz- 
ing manifold shrinks to a point- the number of nearby 
low energy modes is small and our fluctuation integral 
converges in d > 2. This is in accord with the well 
known finite temperature phase transition of the near- 
est neighbor Heisenberg ferromagnet: T c (q = 0) = 0(1) 
(or when dimensions are fully restored- it is of the order 
of the exchange constant). In many rotationally symmet- 
ric incommensurate spin systems (such as the Coulomb 
Frustrated Ferromagnet parameterized by the Coulomb 
strength Q = q A ), a discontinuity in T c may occur as 
Mq — > Mq = q = (q = 0). We anticipate that small 
lattice corrections (A ^ 0) in uq(fc) to yield insignifi- 
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cant modifications to T C (Q): one way to intuit this is 
to estimate T c by the temperature at which the fluctu- 
ations, as computed within the quadratic Hamiltonian 
(6S 2 (x = 0)) =0(1). 

To summarize, we argued that in (essentially hard 
spin) Hcisenbcrg realizations of sufficiently frustrated in- 
commensurate models, no conventional long range or- 
der is possible in the continuum limit at finite tempera- 
tures. More specifically, it is suggested that if the inte- 



gral / 



d a k 



diverges then no long range order is pos- 



v(k) — v(q) 

siblc at finite temperature, unless an expansion about 
the ground state is void (the thermal dependence of 
the higher energy mode occupancies cannot be neglected 
even at T = + ) . If lattice effects are mild then T c (q > 0) 
is expected to be small. In the unfrustrated (q = 0) 
Heisenberg ferromagnet in d = 3: T c = 0(1). Fusing 
these facts, a discontinuity in T c in all rotationally sym- 
metric incommensurate spin systems: 



5T C 



lim[T c (0) 



T c (q)] 



(107) 



is suggested by soft-spin analysis. Such a discontinuity 
will occur as Q — > in the Coulomb frustrated ferromag- 
net. 



XIII. 0(N > 4) FLUCTUATIONS 

We will now find that, on the level of lowest order 
spin wave fluctuation analysis, algebraic long range or- 
der might be possible for n = 4 component spins even 
when their n = 3 counterparts may exhibit exponen- 
tial decay of correlations at all finite temperatures for 
the same Hamiltonian. The fluctuation analysis of any 
0(n > 2) system about a spiral ground state is qual- 
itatively similar to that of the Heisenberg system. As 
noted previously, poly-spiral states will tend to domi- 
nate at large n. The reader can convince him/herself 
that for even n with a p < n/2 poly-spiral ground state 
and for all odd n the fluctuations will give rise to a lead- 
ing order e _1 divergence. The reasoning is simple: the 
poly-spiral states extend along an even number of axis. 
If n is odd then there will be at least one internal spin 
direction i along which S? = and our analysis of the 
Heisenberg model can be reproduced. The lowest eigen- 
energy associated with the fluctuations \ip m ) in the (2p) 
dimensional space spanned by the ground state is higher 
than the lowest eigen-energy for fluctuations along an 
orthogonal direction. The term of constraint is positive 
definite, (ip\H so f t \ip) > 0. If \5S\ <C 1, this implies that 
the quadratic term in 6S(x) stemming from H\ is non- 
negative definite. For Sf(x) — 0, this quadratic term 
in 8Si(x) is zero. The eigenvalue A m ; n = A_ = 2 corre- 
sponds to the zero contribution in 0(5S 2 (x)) from H\. 
Once again, min m {_E m } > Ej &M , from the simple 



mm(^\[H + H soft }\4>) > 

i> 

nmx{<f>\H \4) + mm(^\H soft \0. (108) 

4> £ 

The fluctuations of even component spin about a p = 
n/2 poly-spiral ground state are more complicated. As 
before, coupling between different modes occurs. In this 
case they are more numerous. For n = 4, the fluctuation 
energy, to quadratic order, about a bi-spiral reads 

5H = ^EfcNfc) v(q)]\sm\ 2 + 

^{Efc 1 ,fe 2 [(^l(fcl)^3(fc2) + ^ 1 (fc 2 ) ( 55 3 (fc 1 )) 

X [^<Jl+<?2 + fel+fc2,0 + ^<fl+<f2-fel-fc2,0 + %i+fe2+92-<?l,0 + 
%i+fc 2 +<fl-<j2,o] + 

+ (<ys 2 (fe)tfS4(£i) + 5s 2 (k)ss i (k 2 )) 

X ^k*i+k 2 +qi-q 2 ,0 + %i+£2+g2-9i,0 _ %i +k 2 +qi +<?2 ,0 
fci+fc2— gi— g*2,ol 

- tiSSifySSiih) + 5Si(fci)*S 4 (fe)) 

x [^i+fe+gi+gi.O - %i+£ 2 -gi-g2,0 + ^i+fe+gfe-gi.O 
fci+fe— g2+gi,oJ 

- i(5S 2 (k 2 )5S 3 {k 1 ) + 5S 2 {ki)SS 3 (k2)) 

x fei+fe+gi+gi.O _ ^fei+fe-gi-gi.O + %i+fe 2 +gi -g 2 ,0 
%i+fc2+g2-gi,o]]} 



2ua 
'~N 



L Efe 1 ,fe 2 {^l(fcl)^l(fc2)^ 1+£2i0 + |(% 1+ fe 2+2fl;0 + 



+ 5S 2 (h)8S 2 (k 2 )[5^ +n2fi - Hh 1+ k 2+ 2 ?1 , + 

%i+fa-2gi,o)]} 

+ H #E fc - 1 , fc - 2 {^3(fci)<55 3 (fc 2 )[ ( 5 fc - i+£2 , + i(% 1+ £ 2+29 - 2 ,o + 

%i + fe2-2rf 2 ,o)l 



( 55 4 (/ci)<5S , 4 (fc 2 )[(5j i+ fc 2 



%i+fe 2 -2g 2 ,o)]} 



2 (%i+fe : 



2+2q 2 ,a 



+ 



55i(fci)(55 2 (fc 2 )](% i+ g 2+2fi - % 1+ £ 2 _ 2 ^ ) 

^ ! #Efc 1 ,fe 2 [^3(fc 2 ) ( 55 4 (fc 1 ) + 

5S 3 (ki)SS 4 (k 2 )}(6j: i+j : 2+2?2 - % 1+ £ 2 _ 2 ^ ). 

We may obtain an equation for the eigenvalues by trun- 
cating the expansion for det[H — E] at 0(u 2 ). Just as the 
dispersion relation for low lying states could have been 
derived by setting &ct\H — E] = to 0(u 2 ) and solving 
the resultant quadratic equation for the low lying energy 
states: k = q + S with small S. Here we may do the same. 
The characteristic equation is easily read off: 



N 



= det[W -E} = Y[[(v(ki) - E) 2 } - 

i=i 

4u 2 a 2 Y^[(v(k jl ) + Su)(v(k h + 2qi) + 8u)] 
h 

xl[[(v(h)-E] 2 ] 
-Au 2 a 2 2 J2[(v0 j2 ) + 8u)(v{k J2 + 2q 2 ) + 8u)] 
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x H[(v(ki)-E} 2 } (109) 

where we have shifted £ by a constant. Notice that 
decoupling trivially occurs - terms of the form [v(k) — 
E)][v(k') — E] where the modes k — k' = qi ± q 2 , can- 
cel. The four coupled polarizations break into two pairs 
and that Hi plays the role of a pairing interaction. 
Schematically, for a high dimensional minimizing mani- 
fold M, with low lying states i.e. for the terms containing 
ki = <fi + Si and k 2 = qi + $2'- 

E = E min + al[A n tf\. 1 +A ± 8i. 1 ] 

+a 2 2 [A ll Sl 2 + A ± 5% 2 ] (110) 

with £|| ; _i_ ;TO parallel and perpendicular to the minimizing 
manifold M &t q m , trivially satisfies det[H — E] = 0. 

This dispersion relation agrees, once again, with the 
result derived from the spin wave stiffness analysis 

*E=±J2(v(k)-v(q))\S(k)\ 2 (HI) 
k 

and by expansion of AH for different sorts of twists, the 
dispersion relations of the two spiral simply lumped to- 
gether. When d — 3, as in the 0{2) case (p = 1) this 
dispersion gives rise (in the Gaussian approximation) to 
diverging logarithmic fluctuations: 0(|ln e|). Apply- 
ing equipartition, the Gaussian spin fluctuations in the 
[2i — 1, 2i] plane: 

&Sf 2l _i M (x = 0)> ASf ow [2t _i M (x = 0) 

_ f d 3 k 1 

- kBT J j2nT allA^. + A^i,]- < 112 > 

As the spin fluctuations in the [12] and [34] plane de- 
couple, we may easily compute the correlation functions 
to find that they are a sum of two decoupled pieces of 
the form (Eq.(61)) found earlier for XY systems. For 
general rotationally symmetric incommensurate spin sys- 
tems having a shell of minimizing modes of radius q > 
in momentum space, fluctuations about a certain polyspi- 
ral state explicitly lead to the soft-spin 

4afd 2 cos(gx l; ||) a^Q 1 / 4 
G(x)~ 2^ 2 Li -cxp[-2?77-?7£; 1 (— )], 

i=l,2 -L ZXl 'W 

where ir, n denotes the spatial coordinate in a direction 
parallel to qi etc. As before, r\ = ^^-Q 1 ^ 4 , d — where 
A is the ultra violet momentum cutoff, 7 is Euler's const., 
q = Q 1 / 4 , and E\{z) is the exponential integral. 

For all odd n and for all even n with p < n/2 there will 
be divergent fluctuations similar to those encountered for 
the 0(3) model, no algebraic long range order is found. 

In conclusion, if frustrating interactions cause the 
ground states to be modulated then the associated 



ground state degeneracy (for n > 2) is much larger by 
comparison to the usual ferromagnetic ground states. For 
even n we have found that, generically, the a three dimen- 
sional system will not have long range order when M is 
two dimensional. When n is odd the system will never 
show long range order if M is (d-1) or (d-2) dimensional. 
If, in the continuum limit, the uniform (ferromagnetic) 
state is higher in energy than any other state then, by 
rotational symmetry, the manifold of minimizing modes 
in Fourier space is (d — 1) dimensional. In reality, small 
symmetry breaking terms (e.g. A 7^ in VQ(k)) will al- 
ways be present- these will favor ordering at a discrete 
set of {i^nli^Hi- If n > 2\M\ then, irrespective of the 
even/odd parity of n, then will be a 

( A^- = 0))>,„- 2 ,M!,/^_W_. („,, 

This contribution is monotonically increasing in n; 
Within our scheme, T c is finite and may be estimated 
by the temperature at which the fluctuations are of or- 
der unity. By tweaking the symmetry breaking terms to 
smaller and smaller values, the fluctuation integral be- 
comes larger and larger. For instance, if take A <SC 1 in 
VQ(k) then the integral is very large and T c extremely 
low (in can be made arbitrarily low). Thus as the system 
will be cooled from high temperatures, it might first un- 
dergo a Kosterlitz-Thouless like transition at Tkt to an 
algebraically ordered state and develop true long range 
order at critical temperatures T c < Tkt- 

XIV. THE LARGE N LIMIT- FLUCTUATION 
SPECTRUM, GROUND STATES, GROUND 
STATE ENTROPY, AND T c 

So far we have seen that the classical even n systems 
are more "gapped" than their odd counterparts. This 
is reminiscent of the appearance/non-appearance of a 
gap for quantum spin s chains with even/odd values of 
(2s + 1) [37], [38] and the physics for an even/odd num- 
ber of legs in an s = 1/2 quantum spin ladder [39] where 
such even/odd discrepancies are triggered by Berry phase 
terms. There is no problem in the classical large n (or 
spherical model) limit (just as there is none in the quite 
different quantum spin systems). In this limit, wherein a 
single normalization constraint is imposed 

J2s 2 (x)=N, (114) 

X 

the effective number of spin components n is of the order 
of the number of sites in the system N. The span of the 
system N (the number of Fourier modes allowed within 
the Brillouin zone) is always larger than the number of 
minimizing modes {ql}- In such a case we will be left 
with a divergence as in equation (113) due to the many 
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unpaired spin components. Within the spherical model, 
which is easily solvable the fluctuation integral exactly 
marks the value of the inverse critical temperature 



k B T c 



Ib.2 



d d k 



(27r) d v (k)-v{q)' 



(115) 



Thus, T c = if the latter integral diverges and our circle 
of ideas nicely closes on itself. In terms of ground state 
degeneracies, in the large n limit, any configuration 



S 9 (x) 



E 5 (?> 



(116) 



satisfying the "reality" condition (S(x) — S*(x)) in mo- 
mentum space S*(q) = S(—q), and global normalization 
(Eq.(114)) 



E i^ii 2 



= N 2 



(117) 



is a ground state. The proof is the same as before: any 
configuration S 9 (x) having as its non- vanishing Fourier 
components only those momenta q that minimize the in- 
teraction kernel v(k), is a ground state. Thus, with the 
ground state degeneracy denoted by g, 



ln(5)~|M|. 



(118) 



The ground state entropy is given by the number of min- 
imizing modes in k space- the size of the manifold M 
spanned by these modes. As the spherical model is the 
least restricted of all 0(n) variants, this serves as an up- 
per bound on finite n, 0(n) ground state degeneracies. 
The astute reader will note that the Ising ground state 
constructed in section (IV) for the kernel v z = zA 2 — A 
with z = 1/8 on the square lattice indeed do saturate 
the spherical bound. In the situation that the minimiz- 
ing modes lie on (d—1) dimensional manifold M, the zero 
temperature entropy is sub-extensive yet very large, scal- 
ing with the surface area of the system- \n(g) <~ L d_1 . We 
note that arbitrary periodic patterns (including amongst 
many others Wigner crystals) will be found in the spher- 
ical re-incarnation of many spin models [42] . 



XV. PERMUTATIONAL SYMMETRY AND THE 
INTEGRABILITY OF THE SPHERICAL MODEL 

In many integrable models, the system possesses an 
extremely large symmetry in the thermodynamic limit. 
Perhaps closest in spirit to the spherical model which 
possesses an N\ fold permutational symmetry is a permu- 
tational like character and symmetry within integrable 
Bethe ansatz models where not only the total momentum 
is conserved but also the momentum of each individual 
particle does not change apart from permutations upon 



scattering. Although trivial, the permutational symme- 
try of the large n model has often been overlooked. In 
fact, several authors have attempted (unsuccessfully) to 
find systems having a different geometry of minimizing 
manifolds in k— space (spherical surfaces versus sheets 
etc.) yet with the same degeneracy that have different 
thermodynamics. This quest was unsuccessful for funda- 
mental reasons. The classical spin spherical model (or 
0(n — > oo)) partition function 



const 



n 



i 



(119) 



where the chemical potential \x satisfies 



(2n) d v(k) + » 



is invariant under permutations of {v{k)} 
the above, the permutations 



(120) 



{v(Pk)}. In 



(121) 



correspond to all possible shufflings of the N wave- vectors 
fcj. Although quite simple, this is not universally realized. 
Several authors have attempted to compute the critical 
exponents in the spherical limit (via an RG calculation) 
for systems having different minimizing manifolds yet all 
sharing the same relevant density of states. This quest 
was not very economical. As unrealized by these authors, 
by permutational symmetry these models are identical. 

This simple invariance allows all d-dimensional trans- 
lationally invariant systems to be mapped onto a 1- 
dimensional one. Let us design an effective one dimen- 
sional kernel V e ff(k) by 



/ 



S[v(k) - v]d d k 



,dV f 



eff 



dk 



eff 



(fe)=«- 



(122) 



The last relation secures that the density of states and 
consequently the partition function is preserved. For the 
two-dimensional nearest- neighbor ferromagnet: 



,dV f 



eff , 



= Cl i vT= 



dk 

dx 



1 = P(Veff) 



■ 2 ^i-{v eff + x-2y 



(123) 



and consequently 



\J~Au- 



-)du, 



(124) 



where F(t, s) is an incomplete elliptic integral of the first 
kind. Eq.(124) may be inverted and Fourier transformed 
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to find the effective one dimensional real space kernel 
V e ff(x). We have just mapped the two dimensional near- 
est neighbor ferromagnet onto a one dimensional system. 
In a similar fashion, within the spherical (or equivalcntly 
the 0(n — ► oo)) limit all high dimensional problems may 
be mapped onto a translationally invariant one dimen- 
sional problem, ft follows that the, large n, critical ex- 
ponents of the d dimensional nearest neighbor ferromag- 
net are the same as those of translationally invariant one 
dimensional system with longer range interactions. We 
have just shown that a two dimensional 0(n 3> 1) system 
may has the same thermodynamics as a one dimensional 
system. By permutational symmetry, such a mapping 
may be performed for all systems irrespective of the di- 
mensionality of the lattice or of the nature of the in- 
teraction (so long as it translationally invariant). This 
demonstrates once again that the notion of universality 
(with dependence only on the order parameter symme- 
try, dimensionality etc.) may apply only to the canonical 
interactions. 

Permutational symmetry is broken to 0{f3 A ) for finite 
n. Upon performing a Hubbard-Stratoniwich transfor- 
mation, for a constraining term (e.g. ^Z- In [cosh [77(2;)]] 
for 0(1) spins) symmetric in {i](k)} to a given or- 
der, one may re-arrange the non-constraining term 



En^OPAOMPA;)! 2 ) and rela- 



bel the dummy integration variables H[{rj(k)}] — > 
7J[{?7(P _1 A;)}] to effect the constraining term augmented 
to a shuffled spectra v(Pk). 



XVI. 0(N > 2) WEISS MEAN FIELD THEORY OF 
ANY TRANSLATIONAL INVARIANT THEORY. 



In this section the critical temperature is computed 
exactly for all translationally invariant incommensurate 
0(n > 2) systems. Unlike the situation for Ising spins, 
we proved that simple spiral states are the only ground 
states for incommensurate systems. Here, for 0(n > 2) 
when T <T C : 



(S(x)) = s s 9round - state (x). 
For the particular case 

ground-state^ = ^ . f) 

S° round - state (x)=sm{q-x) 

^ground— state q 



(125) 



(126) 



Now only the ±<f modes have finite weight. Repeating 
the previous steps of section(VI), 



"£V(x = 0,y)S a 2 



round — state 



(y) = 



(127) 



M[Z] EE --H^hr^In^-liz)], (128) 

with [I n /2-i{z)} a Bessel function. With this definition 
in hand, the mean-field equation reads 

|<Si(* = o)>| = s = m[ \J2 v (x>y)( s i(y)) \ ]■ 



The onset of the non-zero solutions is at 



\P c v(q)\ = n. 



(129) 



If V(x = 0) = (no on-site interaction), then 
Jd d k v(k) = 0, implying that v(q) < and T c > 0. 
Note that within the mean field approximation, T c is a 
continuous function of all the parameters in the Hamilto- 
nian. Here the ground state is symmetric with respect to 
all sites. The above is the exact value of T c within Weiss 
mean field theory for the helical ground-states. For poly- 
spirals we will get p identical equations: both sides of the 
self consistency equations are multiplied by af where a/ 
is the amplitude of the I — th spiral in the {{21 — 1), 21] 
plane. As v(q m ) = v(q), we will arrive at the same value 
of T c as for the case of simple spirals. 



XVII. EXTENSIONS TO ARBITRARY TWO 
SPIN INTERACTIONS 

Any real kernel V(x, y) may be symmetrized [V(x, y + 
V(y, x)]/2 — > V(x, y) to a hermitian form. Consequently, 
by a unitary transformation, it will become diagonal. 
The Fourier modes are the eigen-modes of V when it 
is translationally invariant. We may similarly envisage 
extensions to other, arbitrary, V(x, y) diagonalized in an- 
other complete orthogonal basis {|u)}: 



(ui\V\uj) = 6ij(ui\V\ui) 



(130) 



Many of the statements that we have made hitherto 
have a similar flavor in this more general case. 

For instance, the large n fluctuation integrals are of 
the same form 



d d i 



1 



v(u) - v n 



(131) 



with the wave-vector k traded in for u. Here, w m i n = 
min^ v(u). Once again, one may examine the topology 
of the minimizing manifold in u space. If the surface if 
(d — 1) dimensional and v(u) is analytic in its environs 
then, for large n, T c = 0. More generally, in the large n 
limit, the partition function 



Z = const ( TT — ; 



and \(S(x = 0))| = |(5i(x = 0))|ei- We now define 



where the chemical potential \i satisfies 



(132) 
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v(u) + [i 



(133) 



The topology of the ground state sector of 0(n) models 
will once again be governed by a direct product of the 
topology of the minimizing manifold in u space with the 
spherical manifold of the 0(n) group. In the general case 
it will be dramatically rich. 

We may similarly extend the Peierls bounds of 
section(V) to some infinite range interactions also in this 
case by contrasting the energy penalties in the now diag- 
onalizing u basis with those that occur for short range 
systems which are diagonal in Fourier space. For in- 
stance, wc may examine spin configurations computed 
with the non-translationally invariant ( "disordered" ) ker- 
nel V and compare the energies relative to those in a sim- 
ple reference short range interaction which is diagonal in 
k — space (v s hort{k)) and has a minimum v short min = 
minj:{v s i lor t(k)} . The energy penalties for excitations in 
a system with the kernel V, 

AE dlsordered ^ ^^2\S(u)\ 2 [(u\V\u) - v min (u)] 



> 



Vshort 

(k) 

^short min ] (134) 



2N 



if the kernels V and V s hort share the same realizable low- 
est energy Ising eigenstate. 



XVIII. O(N) SPIN DYNAMICS AND 
SIMULATIONS 

We now investigate the dynamics in all (commensurate 
or incommensurate) continuous spin systems. It is seen 
that the equations of motion are relatively trivial in k 
space. Potentially, this allows for the construction of 
new efficient algorithms. If the Hamiltonian 
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then the force, F (z) 



dH 



dS(z 



— , on given spin at size z is 



p & = -\ y) + v & ( 136 ) 



Let an arbitrary constant A satisfy A > — min^{v(k)}. 
The equations of motion in a shifted system with the 
interaction kernel (v(k) + A) are 



d 2 S(k) 
dt 2 



= -[A + v(k)]S{k) 

EE -U 2 J{k), 



(138) 



which may be trivially integrated. Adding the constant 
A merely shifts H — ► H + A/2 with no change in the 
underlying physics. The solution to Eq.(138) is 

-> -> -> dS(k, t) 
S un {k, t) = S(k, 0) cosujkt H ^ — |t=o x uj^ 1 sina^i, 

S U n(k, t + St) = S(k, t) cos u>k6t 

SS(k,t) _! . 
H — — Lo k sincjfeOT. 

This suggests the following simple recursive algorithm: 
(i) At time t, start off with initial values {S(x, t)}; (ii) 
Fourier transform to find {S(k,t)}; (iii) Integrate to find 
the un-normalized {S un (k,t + St)}; (iv) Fourier trans- 
form back to find the un-normalized real-space spins 
{S un (x, t + St)}; (v) Normalize the spins: 



S(x,t + 5t)= S -^,t + St) 
\S un (x,t + St)\ 



(139) 



(vi) Compute {SS(x, t + St)} = { [S(x, t + St)- S{x, t)]}; 

(vii) Fourier transform to find {SS(k,t + St)}; and finally 

(viii) Go back to (ii). Thus far, we neglected thermal 
effects. To take these into account, we could integrate 
these equations with a thermal noise term augmented to 
the restoring force 



= ~v(k)S(k)+F™°%T,t). 



(140) 



Expressed in this format, the execution of this algorithm 
for continuous 0(n > 2) spins seems easier than that for 
a discrete Ising system. Here the equations of motion 
may be integrated to produce arbitrarily small updates 
at all sites. A prioi, this algorithm might be better than a 
real space more brute force approach whereby effectively, 
the torque equations in angular variables (the spins S(x) 
are automatically normalized) are integrated [43] . 



We may symmetrize [V(x,y + V(y, x)]/2 — ► V(x, y), as 
we have done repeatedly without changing H. For a 
more general two spin kernel V(x, y) which is not trans- 
lationally invariant the Fourier space index k should be 
replaced by the more general u. It is readily verified that 



d 2 S(k) 
dt 2 



±v(k)S(k). 



(137) 



XIX. CONCLUSIONS 

In this article, we discussed 0(n) spin models on a 
cubic lattice in mostly translationally invariant systems. 
These include frustrated systems that have been the sub- 
ject of much attention in recent years. Here we outline a 
few of our results. 
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1) Certain frustrated Ising systems were shown to dis- 
play a huge ground state degeneracy. The ground state 
entropy in these systems scaled as the surface area, il- 
lustrating a "holographic" like effect. Computing the 
energies of various contending states is trivial and eas- 
ily allows us, on a non-rigorous footing, to see if lock-in 
effect may arise and what scaling might be anticipated. 

2) The Peierls energy bound may be rigorously ex- 
tended to many Ising systems in which the interaction 
kernel displays minima at a finite number of commen- 
surate wave- vectors. These systems include some with 
arbitrarily long range interactions. This bound can also 
be extended in certain instances to systems displaying no 
translational invariance. 

3) The ground states of all translationally invariant XY 
models were rigorously shown to be spirals (and only spi- 
rals) unless certain special commensurability conditions 
were met. No other ground state are generically possible. 
The trivial uniform and Neel states correspond to com- 
mensurate wave-number ordering (or spirals of infinite or 
two unit length pitch respectively). 

4) The Weiss mean-field equations of the general 0(n) 
systems were looked at. The mean field transition tem- 
perature T MF for all 0(n > 2) models was exactly com- 
puted. 

5) For 0(n > 3) spins, poly-spiral states are the domi- 
nant ones. These correspond to a hybrid of p < Int[n/2] 
(with Int\\ denoting the integer part) spirals in different 
orthogonal planes. Once again, unless commensurability 
conditions amongst the minimizing modes are met, we 
can easily prove that these are the only viable ground 
states. It is also easy to demonstrate that the poly-spiral 
states with the largest viable p (i.e. p = Int[n/2]) are 
statistically preferred. Moreover, they are more stable 
against thermal fluctuations. 

6) We performed a non-rigorous spin stiffness analysis 
to test the stiffness of XY spins under different sorts of 
external twists. We found that certain frustrated systems 
display a smectic like low energy dispersion. 

7) A thermal fluctuation analysis of soft XY spins for 
arbitrary interactions may be easily done. We found that, 
in the general case, a Dirac like equation arises for the 
two component "spinors" . In general, for incommensu- 
rate momenta, the different momenta are coupled. By 
truncating the equations for small u (the soft spin limit), 
we regained the exact same energy dispersion attained by 
the simple minded spin-stiffness analysis. We found that 
in rotationally symmetric incommensurate spin systems, 
the resulting spectra matches that of smectic liquid crys- 
tals. Liquid crystalline properties are natural for many 
of the frustrated systems appearing when competing in- 
teractions are present. 

8) We generalized the Mermin-Wagner-Coleman theo- 
rem to all translationally invariant d < 2 systems with a 
twice diffcrentiable interaction kernel. Here we did not 
only focus on ferromagnetic or antiferromagnetic states. 



All possible orders were excluded. The connection of the 
denominator to Gallilean boosts was noted. At zero tem- 
perature, the denominator of the integrand in the general 
inequality matches that attained by the Dirac like anal- 
ysis. 

9) The resulting generalized Mermin-Wagner-Coleman 
bounds may be examined in high dimensions (d > 2). We 
showed that if a certain integral diverges, the generalized 
Mermin-Wagner-Coleman bounds disallow us to think of 
the dispersion relation for spin waves at arbitrarily low 
temperatures T = + as that computed by assuming a 
perfectly ordered ground state. We also formally linked 
the integrand to a dccohcrcncc (or bandwidth) time scale 
[34] in the quantum case. We further noted that ex- 
actly such a divergence of "decoherence time" is linked 
to a non-rigorous replica derivation of glassy behavior in 
many systems. A large degeneracy and near degener- 
acy may be the underlying cause of two effects- a greater 
fragility of order and sluggish glassy dynamics. 

10) The effect of thermal fluctuation in the most gen- 
eral translationally invariant 0(n = 3) (Hciscnberg) sys- 
tems was considered. We showed that, for incommen- 
surate wave-numbers, the effects of fluctuations here are 
expected to be larger than in XY spins. The only ap- 
proximation made here was that the spin fluctuations 5S 
about the only viable ground state (proved in an ear- 
lier section) were small. This in turn led to a paradox 
at finite temperature: assuming only small spin fluctua- 
tions led to a divergent thermal fluctuation disallowed by 
the normalization of the spin. The parameter u enforc- 
ing the normalization of the spin via the additional term 
Hi = u^ s (S 2 (x) — l) 2 , may be arbitrarily large in this 
treatment (we are not confined to only soft spins 

The only assumption made in this non-rigorous analysis 
that | SS | <C 1. For most canonical systems with kernels 
analytic about point minima in k space (e.g. the nearest 
neighbor Heisenbcrg model), this integral converges and 
estimated T c is sizable. However, for the Coulomb frus- 
trated model and several others, this was seen to suggest 
that T c might be zero as the corresponding fluctuation 
integral diverges in that case. 

11) We extended our thermal fluctuation analysis to 
all translationally invariant 0(n > 4) systems. For even 
n, the determination of the spectrum led to the analysis 
of the Dirac like equation obtained for the XY model. For 
odd n, one component was left unpaired and the fluctua- 
tions were seen to be much larger. The integral appearing 
in the generalized Mermin- Wagner inequality is a strict 
lower bound on the fluctuations. Thus, spin systems 
with a non-trivial minimizing manifold in k space, much 
like spin ladders, display an interesting even-odd effect. 
When we analyze the fluctuations of spins having an even 
number of components (n) , we find precisely the Mermin- 
Wagner integral as the relevant thermal fluctuation inte- 
gral. For a system of spins having an odd number of com- 
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ponents, we find a more divergent (for on-shell minima in 
Fourier space) fluctuation integral. Under the influence 
of the quartic Hi term, spin (or field) components bunch 
up in pairs. For an odd number of spin components n, a 
lone unpaired spin component can give rise to more di- 
vergent fluctuations. Computing the correlators within 
truncated soft spin models, we found that, in rotationally 
symmetric incommensurate systems, weak algebraic long 
range order may persist in even n systems, while being 
absent in odd n systems of the same Hamiltonian. Much 
as in quantum spin ladders and chains, the energy spec- 
trum is more "gapped" and inhibited for even spins. In 
the large n limit, both tend to the "odd" n spherical limit 
fluctuation integral which is none other than the spheri- 
cal limit expression. The bottom line is that the bound 
that we derive by a generalized Mcrmin- Wagner inequal- 
ity is a strict lower bound on the fluctuations. For odd 
component spins n, the thermal fluctuations are much 
larger. This goes against some of the common held intu- 
ition (correct for conventional commensurate ordering at 
q = 0, (n,TT, ...,7r)) that as n is monotonically increased, 
the system orders at lower and lower temperatures as en- 
tropy effects become ever large. By contrast, we find that 
for incommensurate orders, the size of the thermal fluc- 
tuations is not monotonic in n but rather exhibits novel 
even-odd alterations. 

12) The large n limit, which is exactly solvable, gives 
a fluctuation integral (1/T C ) similar to that for the odd 
n case. We further rigorously illustrate how the ground 
state entropy scales as the size of the manifold of minimiz- 
ing modes in k space. In particular, for many frustrated 
models, this leads to an entropy which increases as the 
surface area of the system, leading to a "holographic " like 
effect. 

13) We illustrated how we may think about the topol- 
ogy of the minimizing manifold and use other concepts 
that we introduced hitherto also when the interaction is 
not translationally invariant, as in disordered systems. In 
these cases we considered the new coordinate u param- 
eterizing the eigenstates of the interaction kernel. The 
large n analysis proceeds as before. We illustrated how 
Peierls' bounds may be constructed in certain systems 
with long range interactions. 

14) The dynamics of 0(n) spins in k space was con- 
sidered in the general case, and new algorithm was sug- 
gested. The dynamics in k space is no less easily cap- 
tured than in most standard methods that work only in 
real space. No torques or forces need be considered. The 
evolution of the spin in k space under the influence of the 
two spin interactions is trivial. 
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G(x - y) = (S(x) ■ S(y)) = (cos[0(aQ - 0(y)]>. (141) 

Here 9{x) = q ■ x + A0(af), i.e. A9 denotes the phase 
fluctuations about our spiral ground state and 

G(x-y) =cos(<?- {x-y)){e l{AB ^- AB ^). (142) 

In our harmonic approximation {S6(x)} arc random 
Gaussian variables and only the first term in the cumu- 
lant expansion is non-vanishing. The correlator 



G(f)=cxp[--[(|A0(x)-A0(O)| 2 ])] 



= cxp 



f d d k 1 — cosq-x 
B J (2^F A^ofTAJi 



(143) 



Now let us shift variables k — > k — q = 5, and for purposes 
of convergence explicitly introduce an upper bound on 
< k± < A 



I(x±,x\\) 



1 — cos(q ■ x) 



(144) 



This may be computed by first integrating over fc| em- 
ploying 



f 

j — ( 



1 - cos[a(b - x)] , 7r 

i y n dx = -[1-e ac cos(a6)], (145) 



to obtain 



x z + c z 



A 



II Jo 



1 — cos(fciia:|| + kj_ ■ x±) 
kf. + (A ± ki/A n ) ' 



1 7T / 1 

2 fci A/ Ah^l 



[1 — exp(— \l k\x\\) cos(kj_ ■ x\\ 



(146) 



If <j) denotes the angle between k± and x± then 



2lT 




A± - 

1 — cxp(— < / — — fc_|_ a? 1 1 ) cos(fcx • x±)]d(j) 



k\x\\) / cos(fc^x^ cos (j))dcj). (147) 
Jo 



J {x) 



I 



cos(cccos 



(148) 



the integral of Eq.(144), 



I(x) 



— 4— 7T27T / 

J 



a 1 - exp(-. /j±k' 2 i_x\\)J (k±x±) 



-k±dkj 



Here we follow the beautiful treatment of [44] . Within 
the (hard spin) fully constrained XY model: 



We may now insert the series expansion of Jo (a;) and 
integrate term by term. The Bessel function 
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n=0 y ' 



(149) 



Comparing to the series for the exponential integral 



E 1 {z) = - 1 -\nz-Y J {-T- ( 



n=l 



n(n!) 



(150) 



(7 is Euler's constant), we find that 
G(x) ~ — exp[— 2i]-f — r]Ei(- 



= )] x cos[qx||], (151) 



4 n 

where 77 = = if§-? (in the last equality Eq.(60) 

is invoked) , and d = =^ , with A the ultra violet momen- 
tum cutoff. 
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More generally, in d > 2 dimensions, the leading order 
infrared contribution reads 
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the average charge density over the entire system (the 
background charge taken to be a half in the standard 
Ising lattice gas problem). Within the spherical model, 
only global normalization of the spins is required and the 
charge density may be arbitrary at any given site. As the 
Hamiltonian is quadratic in {Si}, the average (S) = 0, 
self-consistently implying that the average charge density 
(p) = P- 

[43] In real space, the equations of motion trivially read 



d 2 



dt 2 



^V(£-y)sin[#H)-0(y)]. (159) 



for an 0(2) system. For the three-component spin system, 
with 9 and 4> the longitudinal and azimuthal angles, 

— V(x — y) s'm9(x) x s'm9(y) s'm[4>(x) — 4>(y)]; 

v 

d 2 e(x) i . d^g) . 2 

= -sm29 (x)( — — L ) 

dt 2 2 y Jy dt ' 

+ ^Tv(x-y){sm6(x) 

y 

cos 6{y) - cos 6»(f) sin 9{y) cos^(f) - (p(y)]}. (160) 
[44] J. Als- Nielsen et al., Phys. Rev. B 22, 312 (1980) 
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